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ABSTRACT 

The monitoring of the deformation response of the rock 
mass upon excavation of a tunnel is the most commonly 
applied means of evaluating the design and stability of the 
Opening. Because of the importance attached to this 
monitoring operation, a critical assessment of current 
practice, with particular regard to excavations in highly 
stressed or overstressed rock, was undertaken. 

Laboratory model tests were employed to determine the 
rock mass response to the excavation of a tunnel under 
controlled environmental conditions. These tests were 
classified as either excavation simulation or external 
loading. In the former, a tunnel is excavated or enlarged in 
a pre-stressed sample while in the latter, the sample 
containing a pre-excavated tunnel is loaded at its 
boundaries. Plane strain conditions were maintained along 
the the tunnel axis in all tests. In excavation simulation 
testing, the development of radial strain around the opening 
was recorded during and following excavation. The 
convergence of the tunnel walls could only be monitored in 
the post-excavation period. In external loading tests, both 
convergence and radial strains were measured throughout. 
Coal was chosen as the model material because of its small 
scale, regular joint pattern, low strength and 
time-dependent deformation properties. 

As a result of the model tests, a rational procedure 


for the evaluation of tunnel behaviour, based on the rate of 
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tunnel convergence, was developed. The appropriate rate 
criteria for this evaluation require field determination and 
their application requires separation of the convergence due 
to face advance from that due to the creep of the rock mass. 
Expressions for the determination of the convergence and 
radial strain rates due to the advancing face were developed 
for this purpose. 

The limitations of current deformation monitoring 
practice in tunnels constructed at great depth were 
identified. It was determined that the interpretation of the 
results can be signifigantly influenced by the location and 
time of installation of the instrumentation, particularly 
where the material exhibits time-dependent behaviour. To 
evaluate the design and stability of the excavation requires 
multi-directional instrumentation installed as close to the 
face aS practically possible (preferably ahead of the face, 
although this is generally not possible at depth) with 
frequent reading and consideration of the face advance 
effect. It was further determined that there is a need for 
instrumentation capable of displacement detection ahead of 


the face. 
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1. INTRODUCTION 

The greatest demands of modern societies are the same 
as those of ancient societies: resources and accomodation. 
While past generations were able to avail themselves ina 
relatively easy manner (though not necessarily easy given 
the technological level) the increased demands of a growing 
world have caused 'modern' men to foray in further fields. 
Past exploitation of surficial resources have forced man 
deeper beneath the earth's surface for the required 
commodity. Land and energy costs have halted the post-war 
urban sprawl and forced man into vertical accommodation, his 
living space upwards and storage space downwards. These 
resource and storage requirements, along with their support 
measures are the catalyst for the field of underground 
excavation (tunnelling) technology. 

While underground excavations have been constructed 
throughout recorded history, it is only recently that both 
the scale and economics of tunnelling have undergone 
dramatic change. Prior to this current, cost-conscious 
efficiency, tunnelling was considered an ‘'art' practiced by 
a few 'artists' who could cope with the existing demand. 
Present demand, however, requires that the behavior of 
underground excavations be understood so that less 
experienced engineers or contractors may undertake these 
endeavours, in essence the transformation from an art to a 


science. 
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Current tunnelling technology, while far advanced, is 
Still in this transition stage. Because of this, underground 
excavations, whether they be in soil or rock, and 
independent of their intended use, require monitoring. 
Unfortunately, field monitoring in itself is not sufficient 
for the development of a predictive capability. In the 
field, there are a multitude of influences on tunnel 
performance which complicate analysis. Controlled 
environmental conditions are essential for a rational 
analysis. One means by which this may be achieved is the 
Simulation of field conditions in laboratory scale models. 
In this manner, the response of an excavation to 
pre-determined stimuli may be observed and a clearer picture 
of the overall behavior of underground excavations emerge, 


especially when coupled with field observations. 


1.1 Thesis Objective 

The monitoring of underground excavations 1s crucial to 
the further development of tunnelling technology. Because of 
this importance, it is essential that a proper framework for 
data interpretation be devised, the limitations to the 
instrumentation employed for monitoring be understood and 
that their application be appropriate for their purpose. 

Most tunnel excavation monitoring schemes employ 
deformation instrumentation as the primary mode of 
investigation. This is particularly true in excavations 


where large deformations are anticipated, such as in weak 
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rock or rock exhibiting time-dependent deformation. It is 
therefore essential that the use of current deformation 
instrumentation in this type of environment be critically 
reviewed. That is the objective of this thesis. 

In order that this objective be realized, it was 
decided to forego field analysis, due to its inherent 
ambiguities, and instead, concentrate on laboratory model 
Studies. Particular attention was paid to the evaluation of 
opening stability by deformation measurements and the 


instrumentation response to tunnelling activity. 


1.2 Scope of the Thesis 

The first step in any investigative research is the 
review of the current knowledge on the subject. The results 
of a review of current tunnel monitoring as practiced by 
both the civil engineering and mining engineering 
(particularly coal mining) professions are presentd in 
Chapter 2. Areas of particular interest included the 
monitoring objectives, the means by which these objectives 
are achieved, the selection of the type of monitoring, the 
type of instrumentation, the selection of appropriate design 
and performance criteria and the limitations as presently 
understood. 

Decriptions of the test apparati employed in this 
research are presented in Chapter 3. Also included in this 
chapter are physical descriptions of each sample along with 


a description of the type, installation and layout of 
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instrumentation and a summary of the adopted loading 
sequence. Coal was used as a test material because of its 
time-dependent strength and deformation properties. 

In order to analyze the observed behavior of the 
tunnels during process simulation testing it was necessary 
to determine the material properties of the test specimen. 
Both the time-dependent and time-independent properties were 
obtained from the response of large intact blocks 
(600x600x200mm) to a variety of loading conditions. The 
results of these investigations, presented in Chapter 4, 
allowed formulation of numerical expressions for the 
prediction of the deformation response of the test material 
under given stress conditions. The analysis of the test 
results, with particular emphasis on the implications 
towards deformation monitoring in underground excavations in 
weak rock subject to time-dependent deformation, are 
presented in Chapter 5. Attention was focused on the 
limitations of deformation monitoring for stability 
evaluation and the verification of design assumptions under 
both advancing tunnel and long-term conditions. 

Conclusions regarding the application of deformation 
monitoring instrumentation in weak rock at depth are 
presented in Chapter 6. Also included are recommendations 
for the implementation of deformation instrumentation and 


suggestions for further research in this area. 


+ ee So, 

AN ee coe 7 a Lede wat 0 nots 

ee eee | Haaser ee s 

#3 en oe yom 
| pte mae ike _ a 


_ wit To deine 

yiseesasn 2ew DL esa 

.oamizaqe 2264 si? aoe }12g0 

ssw @eitre@esg sichaaeeninelaot 
eAvoid sasint- spasl h a ao: 

eit. .enerzibmes. pittbaot 20 geeksatealh 


[. 


Lb <#3qsd0 ni bedneesig) enob dene yaapa? i = sity Be eels 7 
ad3 102 anoleasigze Les yssmun” to cosnacumnet Sovota > 
istivetam i423 385 To iat ha Ao etamrivtad ode 20.7 ; _ 


1282 sd Jo efaylons en -2indiaiBdos gait vevig sha 


enotssabiemi sdf ao eteurigns: astucisveq dtiw. 

ni aactievssxs a a ib: onisostnem noiseavanab , 
ers «nots Biol ‘wacbicent-haity ‘em’ d5aheam Me aor 
aia no beeuse? meytnotahet2s) V8. ‘yeeqed od beenseexq, | 
vat i Lidese sm? heal aaa nosdamnbtab io eee 
1sbay ahoi iquess: dob? 2b) oatdeataiagy, ada. Sap f | 
aenanrar™ Pane sane 
aoisenroreb #0) ne 

| xe Aagaty tacdsox nse | 

ariei nena ~ aoe mer resaedo 


— 


\9 
= 


2. REVIEW OF CURRENT DEFORMATION MONITORING PRACTICE 


Zein crodquction 

The 'art' of tunnelling is becoming more of a science 
as a result of the application of recent advances in design, 
construction and monitoring of underground openings. These 
three components are no longer discrete elements of the 
tunnelling process but are interrelated as shown in 
Figure 2.1. In the site characterization phase, the geology 
along the alignment is determined as are the engineering 
properties of the rock mass required for its description. 
Once these are known, the excavation sequence, construction 
method and support system are often selected on the basis of 
Simple empirical classification systems. Detailed design is 
then derived from empirical methods, more refined 
classification systems or an analytical approach. After 
commencement of construction, monitoring of tunnel 
performance permits critical evaluation and possible 
revision of the design or construction procedure. Because of 
the importance of this final phase to modern tunnelling 
technology, it was decided to critically review the means by 
which the performance is presently evaluated. 

The above tunnelling process is generally limited to 
civil engineering works. In mining, the most common design 
concept is the duplication of previous support and geometry 
from existing operations under similar conditions 


[Dahl1(1978)]. Monitoring of underground excavations, 
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Figure 2.1 Tunnel Design Sequence 
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however, has been actively pursued by the mining industry 
and valuable contributions in this field have resulted from 
their efforts. While the purposes of the underground 
excavations for civil and mining works may be very 
dziferentyathe rockhmassoisaunaware of it;sandowill reactoin 
a Similar manner to similar stimuli. The advances in the 
monitoring of underground excavations in rock are equally 


applicable to both. 


2.2 Role of Monitoring in Underground Excavations in Rock 
The monitoring of an underground excavation in rock is 
Simply the measurement of the rock mass response to the 
changes imposed. These measurements may then be employed to 
Satisfy the objectives of the monitoring program. Paramount 
amongst these objectives are the determination of the 
opening stability and the verification of the design [Kaiser 
(1981d),Franklin (1977)]. Additional objectives have been 
defined by Lane (1977), Londe (1977) and Dunnicliff et al 
(1981). These are outlined below along with those previously 


mentioned. 


Stability Evaluation 

The primary role of monitoring is to ensure the safety 
of the underground environment during both construction and 
operation. Current monitoring technology, if applied 
correctly, permits determination of incipient tunnel 


instability with sufficient time for excavation or remedial 
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action if the process leading to instability is gradual. 


Design Validation 

The design of the tunnel support is made with gross 
assumptions as to the geology and rock mass properties. To 
ensure that the design iS appropriate requires verification 
of the assumed conditions and predicted performance by 
in-Situ meaSurement. 

The validity of a new design approach may also require 
the performance monitoring of a prototype for comparison 
with the new and proven design methods. Advancement of the 
State-of-the-art requires hard evidence which can only be 


achieved by in-situ meaSurements. 


Construction Assessment 

The adequacy and suitability of a given or new 
construction method is best assessed by monitoring the 
performance of the excavation during these endeavours and 
then comparing with measurements made in excavations with 
Similar geologic conditions but where different construction 
methods were employed. 

As tunnels become more expensive to construct, the 
avoidance of overdesign is of increaSing importance. In situ 
instrumentation allows the designer to implement a ‘minimal 
support' design with provision for additional support where 


measurements dictate. 
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Problem Diagnosis 

The remedy of an incipient instability requires a prior 
knowledge of the nature of the instability (i.e., the 
mechanism, location and extent). This can only be achieved 


by insitu measurements. 


Performance Documentat ion 

Records of the performance of the underground 
excavation provide documentation of ‘as built' conditions 
and structural response. These may then be employed to check 
for compliance with contract requirements and to assist in 


the design of future works. 


Public Reassurance 

Use of current tunnel monitoring technology serves to 
reassure the workers and the public in general of the 
owner's and contractor's concern with their safety as well 
as with the long term integrity of the excavation. 

The objectives pursued by the monitoring program at any 
given site will depend on the nature and difficulty of the 
project, the owner's habits, the designer's style and the 
contractor's experience [Londe (1977)]. The means by which 
the monitoring is performed are discussed in the following 


section. 
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2.3 Deformation Monitoring 

In addition to visual observation of the excavation, 
monitoring is performed by in-situ measurement of stress 
changes and displacements. The environment in which these 
measurements are made is particularly harsh. Therefore, it 
1s necessary that any instrumentation be of simple and 
rugged design yet yielding reliability and accuracy. 
Readings must be easily interpreted so that the necessary 
action may be taken and neither the instrument nor its 
reading should duly interfere in the construction or 
operation of the opening. Given these operating criteria, 
the most commonly employed monitoring system is that of 
Gisplacement measurement. Londe (1977) considered stress 
change measurements the most difficult to interpret as the 
instrument itself alters the local stress field and the 
discontinuities inherent in the rock mass cause an 
unacceptable degree of scatter in results. Kaiser (1981d) 
considered displacement meaSurements as the most important 
mndicator un) Ground control aS did Cording et aj] (1975) 
because of the ability of deformation measurements to yield 
information as to the depth and magnitude of the movement. 
Lemcoe et a] (1980) in their state-of-the-art review of 
measurement techniques found displacement measurements more 
Girect, easier and less costly to perform and hence more 
commonly employed than stress measurements. The 
"observational method' of tunnel design as exemplified by 


the New Austrain Tunneling Method (NATM) depends heavily on 
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interactive deformation instrumentation [Kuesel (1983), 
Amberg (1983)]. In the mining environment deformation 
measurements are employed in the design of roof and face 
support sfandiuootfortilooraboltings ituy 61982) ) cotairn 
assessing stability [Whittaker (1974)]. 

Due to the dominance of displacement or deformation 
measurements in the performance monitoring of underground 
excavations it was decided to forego any detailed 
conSideration of stress measurements in monitoring. This 
does not imply that stress measurements are not of value in 
monitoring. Indeed, stress change meaSurements can greatly 
enhance the value of deformation measurements [Kaiser and 
Mackay (1982)]. 

The monitoring of deformations in underground 
excavations is achieved by optical surveying, convergence 
measurements and borehole extensometer measurements. 
Generally no single means of monitoring is employed on a 


given project as the three methods complement each other. 


2.3.1 Optical Surveying 

Optical surveying methods such as levelling and 
triangulation are widely employed to monitor movements of 
the rock surface or tunnel lining [Londe (1977)]. Where 
accessible and where the meaSurements can be related to a 
remote stable base the absolute displacements of the tunnel 
wall can be determined. The advantages of this system are 


that the personnel and equipment are generally on site (for 
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alignment control) however, these are outweighed by the 
disadvantages of interference with the construction 
procedure and the considerable time required for computation 
and interpretation. This latter point is particulary 
important where monitoring is employed to detect incipient 
instability. 

Electro-optical distance measuring instruments 
employing a modulated light or laser beam projected onto 
reflective targets on the surface of the excavation are 
Superceding the optical instruments [Franklin (1977), Hoek 
and Brown (1980)]. Their principal advantage over the 
optical system is that many measurements, including those in 
highly inaccessible regions can be made from a single 
instrument location in a short period of time. 

Elaboration on this form of instrumentation will not be 
undertaken as laboratory simulation is not possible. 
However, the results of the laboratory investigation 
reported in Chapter 5 have significance with regard to this 


form of monitoring as well. 


2.3.2 Convergence Measurement 

One of the most common and easily performed deformation 
monitoring methods undertaken in underground excavations is 
the measurement of the change in diameter or convergence 
between reference points. Convergence is measured by either 
a tape or rod extensometer with the tape being the more 


widely used because of its range (50 m versus 8 m for the 
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rod type) and handling ease [Lemcoe et a] (1980)]. 

The tape extensometer consists of a stainless steel 
tape with a tensioning spring and a mechanical measuring 
device such as a dial indicator gauge. The rod extensometer 
employs a rod in the place of the tape, thereby removing the 
need for a tensioning device. The rod can be of fixed length 
or telescopic (possibly spring loaded) and is fitted with a 
micrometer or dial gauge [Franklin (1977)]. Lu (1982) 
describes an automated rod type convergence meter where, 
upon convergence, a lever arm on the device causes a scribe 
point to trace a path on a pressure sensitive chart or 
rotating drum giving a complete time-convergence record. Rod 
type extensometers employing electric readout devices and 
their application to mining are described by McVey and Howie 
(198 19)% 

Convergence targets for both systems generally consist 
of steel or aluminum anchor bolts secured in the rock by an 
expanding shell anchor or by grounting. The target ends are 
designed to mate with the appropriate convergence meter 
[Franklin (1977)]. 

The accuracy of these instruments is highly dependent 
on the care employed in making the reading (temperature 
compensation, spring tension, etc.) and in installing and 
providing protection for the reference target [Lemcoe et al 


(1980)]. 
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2.3.3 Borehole Extensometers 

The previously described deformation measurement 
systems give only surface displacements of the opening 
walls. For 'in-depth' measurement of the rock mass 
Geformation, borehole extensometers are required. As the 
name implies this type of instrumentation requires the 
drilling of boreholes either from the excavation or from 
outside the excavation if conditions are favourable (shallow 
cover or parallel adit). Two types of borehole 
extensometers, probe and fixed, are commonly employed 
[Lemcoe et a] (1980)]. In both systems the differential 
movement of anchor points installed at various depths is 
measured. 
Probe Type 

The probe type extensometer system meaSures the change 
in location along the borehole axis of pre-set anchors with 
respect to each other or some datum (furthest anchor or 
borehole collar) by means of a moving sensor. The most 
popular of this type employs magnetized anchors and a probe 
capable of sensing the magnetic field [Lemcoe et a] (1980)]. 
The major disadvantage of this system is that it is labour 
and time intensive and not suitable for automated or 
continuous recording. 
Fixed Type 

For increased precision or where continuous or remote 
readout is required extensometers of the fixed type are 


used. The displacement along the borehole axis of selected 
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points is determined by measuring the displacement at the 
borehole collar of rods or tensioned wires extending from 
anchors located at those points [Franklin (1977)]. 
Measurements are commonly made with a dial indicator gauge, 
however for continuous or remote monitoring electric 
displacement transducers are employed [Amberg (1983), Lu 
(1982)]. Rod type extensometers yield more reliable and 
accurate results as the wire type are prone to kinking and 
creep. In single hole, multi-point installations, wire 
extensometers are more often used because of their reduced 
Space requirements [Franklin (1977)]. 

Because of the ability of the borehole extensometer to 
provide information regarding the extent and depth of rock 
movement they are generally used in conjunction with 
convergence meaSurements to monitor underground excavations 
during the construction and operation phases. Selection of 
the monitoring instrumentation is only the first phase of 
the monitoring process. When and where to monitor and what 


criteria to employ are of equal if not more importance. 


2.4 Instrument Location and Monitoring Schedule 

The general principle governing the selection, location 
and recording schedule of deformation meaSuring instruments 
is that of working from the general to the particular 
[Franklin (1977)]. The intended role of the monitoring 
system will also influence the selection as will the type of 


excavation and support method [Cording et a] (1975)]. 
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Franklin (1977) suggests that when monitoring is 

- performed to ensure the safety of the excavation, there 

should be extensive coverage with low precision, low 

frequency monitoring supplemented with more frequent, more 
accurate monitoring where conditions dictate. If design 
verification is the objective, he suggests selecting 

"representative profiles’ corresponding to various geologies 

and geometries for monitoring. Borehole extensometers are 

required in addition to convergence measuring devices at 
these profile locations. Other design principles presented 
by Franklin include: 

1. Adoption of redundancy in the instrumentation to 
accomodate instances of damage to an installation or 
malfunction of an instrument; 

2. Employment of a large number of low precision 
instruments rather than fewer, high precision 
instruments. With more instrumentation, anomalies in 
readings are more readily ascertained; 

3. Adjustment of the monitoring frequency to suit the 
Situation. After installation monitoring should be 
frequent to identify trends in ground behavior. 
Afterwards the frequency of readings will depend on the 
activity of the specific site; and 

4, Proper installation and calibration of the instrument so 
that it functions as intended. 

Both Franklin (1977) and John (1977) recommend convergence 


measurements spaced at 50 to 100 metre intervals along the 
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tunnel axis and at geological changes under good rock 
conditions. John also suggessts a reduction in this interval 
to between 10 and 20 metres when fracture processes are 
occurring. At each measuring station the convergence is 
measured at the springline and from the crown to the invert. 
Supplementary measurements across oblique diameters are 
performed when considered necessary. At selected monitoring 
Stations borehole extensometers are also installed, 
principally to verify support design. Cording et al] (1975) 
recommend installation as close to the face as possible so 
measurements can be made before tunnel advance. They further 
recommend that the anchors be placed at a minimum of two 
depths; one close to the tunnel wall (1-2 m) to record near 
Surface movement and the other deeper to give evidence of 
deep seated movements. Extensometers are generally situated 
in the crown and at the springline with occasional 
installation in the invert and at oblique angles. A typical 
deformation instrumenation array is shown in Figure 2.2. 
This array demonstrates the basic needs of the 
instrumentation system for a full face, circular excavation. 
Horizontal and vertical convergence determine the deformed 
shape of the excavation. Extensometers in the crown and 
invert as well as convergence meaSurements from the 
springline to the crown allow separation of the convergence 
components (i.e., bed separation, floor heave, etc.). 
Ideally all deformation measurements are tied together and 


then referenced to a remote datum (by tunnel survey) in 
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Figure 2.2 Deformation Monitoring Instrumentation Layout 
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order to fix the displacement field in space. In multiple 
heading excavations the same principles are employed, 
however, the location of the instruments will depend on the 
particular excavation sequence and the anticipated behaviour 
as discussed by Cording et aj (1975). 

Deformation monitoring in mining applications is not as 
rigorously applied as in civil engineering projects due 
mainly to the short life span of the majority of mine works. 
Where monitoring is employed it generally consists of roof 
to floor convergence measurements. Occasionally the 
convergence is further defined by separation into floor 
heave and roof sag [Hargraves and Martin (1977)]. Stability 
monitoring by microseismic methods has been undertaken by 
the mining industry [Blake (1982)] but is not considered in 


this research. 


2.5 Stability Criteria 

When deformation monitoring is being employed to 
ascertain the stability of an underground excavation, the 
recorded measurements must be compared to design criteria in 
order that an assessment be made. Both Franklin (1977) and 
Londe (1977) consider absolute values of displacement as 
insufficient for stability assessment. They along with John 
(1977) and Cording et al (1975), however, consider the use 
of displacement rate as a sufficient, suitable criterion. 
Indeed this is the most frequently employed means of 


Stability assessment in tunnelling [Lane (1977)]. 
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In mining, closure rate measurements are commonly 
employed in the determination of opening stability [McVey 
and Howie (1981), Hams (1978) and Wieselmann (1968)]. Lu 
(1982) and Pothini and von Schonfeldt (1978) recommend both 
the cummulative convergence and the rate of convergence be 
used in establishing stability. Pothini and von Schonfeldt 
also consider the change in the convergence rate in their 
assessment of opening stability. Serata et aj (1976) suggest 
that the determination of opening stability can best be 
determined by plotting the closure rate versus the 
excavation age in a double log diagram. The degree of 
current stability is indicated by the magnitude of the creep 
rate with that of future stability by the slope of the creep 
rate curve. This point will be discussed later. 

Selection of an appropriate critical deformation rate 
for stabilitycassessment is not “ancveasy task. Thesrange of 
convergence rates encountered in underground excavations is 
quite large. John (1977) found that rates less than 10 mm 
per month (in an 11m dia. tunnel) resulted in low lining 
loads while McVey and Howie (1981) found rates of 5 mm per 
minute to indicate impending roof fall in a specific coal 
mine. Criteria such as these indicate the range of rates 
which may be encountered. They are not directly transferable 
to other sites. 

As an early indication of instability, in the absence 
of predetermined values, high rates of movement unrelated to 


excavation advance or sudden rate increases may be used 
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[Cording et al (1975)]. To determine the site specific 
critical rate requires observation of rates, under similar 
conditions, where stability and instability of the 
excavation were encountered. These can be determined from 
other projects but care must be exercised in any 
extrapolations needed due to changes in support or geometry. 
The best method of selection is the observation of 
deformation rates in pilot headings or test sections of the 
Same project. 

In mining, special care is required to separate the 
effects of headings at other levels or in parallel advances. 
The excavation sequence (i.e., leading, trailing or 


Simultaneous gateroad excavation) is also relevant. 


2.6 Limitations of Current Deformation Monitoring Practice 

Current deformation monitoring practice in underground 
excavations suffers limitations in both concept and 
application. The former relates to the inadequacy of the 
program as discussed by Kaiser (1981d) while the latter 
relates to field problems presented by Dunnicliff et al 
(1981) and others. 

Kaiser (1981d) points out that most monitoring systems 
concentrate on the determination of the stress or 
displacement field while neglecting or inadequately handling 
the determination of the material property distribution and 
the mode of behaviour of the system. Based on the results of 


model studies, Kaiser concludes that interpretation of field 
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behaviour is not possible without a knowledge of all three. 
While acknowledging that it is often necessary and indeed 
permissible to make assumptions for some of the components 
described above, he suggests that this is not the case with 
the mode of behaviour. The monitoring system should be 
designed to ascertain the deformation process rather than to 
record the response to an unknown process. 

Applications problems related to deformation monitoring 
are of both procedural and mechanical origin. Procedural 
problems identified by Dunnicliff et a] (1981) include: 

1. Failure to identify the need for instrumentation at a 
given location; 
2RSUinilexibiiatyiinethe instrumentaion layout and recording 
schedule (ie. failure to adapt to changing 
requirements); and 
Seer Oi latoryeco. lecticnaand interpretation of measurements. 
An additional procedural problem relates to the reduced 
deformations measured by instrumentation installed behind 
the face as noted by Lo and Lukajic (1984). Mechanical 
problems are generally related to instrument reliablity 
particularly where anomolous measurements are recorded. In 
Situations such as this the experience of the 
instrumentation engineer is usually relied upon [Cording et 
al (1975)]. Minor mechanical problems occur due to improper 
installation or post installation damage. 
Additional concerns regarding the monitoring of 


underground excavations were identified by Wieselmann (1968) 
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and Muir Wood (1979). Wieselmann in an investigation into 
the use of closure rates aS a monitor in potash mines found 
that the rates were dependent on the opening shape and the 
method of excavation. Muir Wood indicated that, in materials 
exhibiting time dependent behavior, separation of the 
deformation into that resulting from the variation of the 
Stress field due to face advance and the creep under 
constant shear stress is necessary. 

In an attempt to assess the implications of these 
influences on the deformation monitoring of underground 
excavations, particulary those in weak rock, a series of 
Simulation tests were performed on laboratory specimens. The 


results of this investigation are presented in Chapter 5. 
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3. SAMPLE DESCRIPTION AND TEST METHODOLOGY 


3.1 Sample Description 


Sells Description of Coal 

The coal employed as the model material in this testing 
program was obtained from the Highvale Mine in the Wabamun 
Lake District near Edmonton, Alberta. This coal was graded 
by Standsfield and Lang (1944) as subbituminous B according 
to the Canadian Classification and its composition 
determined by Pearson (1959) as presented in Table 3.1. A 
detailed structural survey conducted by Noonan (1972) 
concluded that the coal is horizontally bedded with two 
orthogonal cleat systems perpendicular to bedding. The major 
System is oriented N45°E and has an average joint spacing of 
approximately 20 mm. Jointing is essentially planar but 


discontinuous. 


3.1.2 Sample Preparation 

Test specimens of approximately 600 x 600 x 200 mm were 
prepared from large block sample as described by 
Kaiser (1979). The major joint set was oriented parallel to 
the main diagonal of the specimen and the bedding planes 
parallel to the larger surface. In this manner, the bedding 
plane was perpendicular to the tunnel axis and the major 
jointing at approximately 45° to the principal stress 


direction when installed in the test frame. 
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Table 3.1 Composition of Coal 


Moisture 

Ash 

Volatile Matter 
Fixed Carbon 
Fuel Ratio 


Gross Ku per N 


18.4 

Oa) 
24.4 
38.9 
50 
1898 


*% by weight 


2116 
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3.1.3 Description of Individual Specimens 


Sig a Ala 1 EET enue Y (eas 

Sample MC-5 was prepared in the prescribed manner to 
final dimensions of 620 x 615 x 210 mm. The major joint set 
was inclined at approximately 49° to the major principal 
stress direction. Joint spacing varied from approximately 10 
to 20 mm. A plan view of the top surface of this sample is 
shown in Figure 3.1. 

A fracture running linearly across the sample from 
Side 2 to Side 4 was noted during pre-test mapping, however 
it was located outside the tunnel area and thought to be of 
minimal significance. Associated with this fracture was a 
large depression caused by removal of loose material. The 
depression was filled with high strength grout which may 


havelmntluenced themourueness ins this region. 


3. i.3'92>Sanple* MeZo 

Sample MC-6 was prepared in the same manner as MC-5 to 
final dimensions of 620 x 620 x 206 mm. The major joint set 
was®d bstinctly_non planar with inclinations “of.21.° to.68° 
from the major principal stress direction. Joint Spacing was 
generally in the range of 10 to 20 mm with spacing of less 
than 5 mm in the heavily jointed region surrounding a clay 
shale inclusion. Plan views of the top and side surfaces are 
presented in Figure 3.2. 

A clay inclusion, elliptic in section, is evident in 


the plan view of Side 3. This inclusion dipped at 


- 
44 
. 


em “sit 03> °@5 ae sone $s bert Ls 
Ot ylossiizotggs moat bedsaveD oni sede sakes prone 22 
ai sidaise efdy 30 Soetdue coe sit) ia vole nag a mnt 08 1 o2 
| ane smett ni e bs 
mort siqmec sdz 226 is vixeants weortignife gare a 

yevewod ,ehagqgqsn 2564q) pn das exer bios abi@ of 

10 4d of sipudrte Gis face fears site etsy 3) maee ia . 
6 aBw 931 Eye aids: dtiw bevei nomae saangzikinpte Semi 


na 


sat beteaaaee se rool io Lavousn “cd aeeuint noigesigeb a ft 
sf; 
yen doidwW Sbo%p dipnsive deid dice atti 2nw notees qob 


.ngige: aang ni neha ide, pity geomaulind eved 7 


. ) + 

hes | plete | 

o2 a-OM 28 weaniem ems2 ed3 ni. 4s ae 25¥ 3-2M ole 
ts2 Infef Bh ait at 8 ‘BOS Os@ x OSby/t0 saa c tenit 
"80 o2 “4S io any OF Are 3 he ‘snustg-hon indie = hae 
26W bn? sa fatot nokioeyis seer te Lwakoning pay ws mor} 
neat 19 pnisbye git mH OS oF Of es) — i a 
vats Bi. encenuore tess bpsaret eh 
8 29382 + 30 vada ut 
2 7 ~ se som bes 
fib on oni 


uy oni ait ae nm ¢ ee oe f 
a ir J = 


> » 


(iis 
as 
m ¢ | 


ve ie 


_ 


we reve Base i 
ff bie bn re | 


7 


of Sample MC-5 


of the Top Surface 


3.1 Plan View 


Figure 
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Figure 3.2 Plan Views of the Top and Side Surfaces of 


Sample MC-6 
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approximately 7.5° towards Side 1 and increased in width 
down dip. Jointing, which in the remainder of the sample was 
essenvularly vertical, “inclined wminethesvicinity Obmthe 
inclusion; the joints radiated from the inclusion. The 
inclusion had a dominant effect on the structure of the 
sample and will thus have a considerable influence on the 
Sample performance. Side 2 and, to a lesser extent, Sides 1 
and: 3, nequiredsashighestrength? grout infilling) of 
depressions caused by removal of loosened material. As these 
infillings were concentrated around the sample periphery, 
their effect on near tunnel performance should be 


negligible. 


Sl. Ses Sample MC-7 

Sample MC-7 was prepared in the required manner to 
final dimensions of 625 x 625 x 212 mm. The major joint set 
was non-planar and non-uniform. The strike of this joint 
set, referenced to the major principal stress direction, 
varied from 20° to 70° with a mean value of approximately 
40°. The dip angle varied from 70° to 90°. The joint spacing 
varied from 5 to 25 mm. Plan views of the top and side 


surfaces of this sample are presented in Figure 3.3. 


3.2 Deformation Instrumentation 

During the testing of a given specimen the following 
deformation measurements were recorded: 
1. boundary displacements; 
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Figure 3.3 Plan Views of the Top and Side Surfaces of 
Sample MC-7 
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3. tunnel closure. 

The means by which these measurements were obtained are 
briefly described below. More detailed descriptions of the 
deformation monitoring instruments are presented by Kaiser 


(1979) and Guenot (1979). 


3.2.1 Boundary Displacements 

The boundary displacements in the principal stress 
directions were measured by displacements transducers 
(LVDTs) mounted on metal stands attached to the base plate 
of the specimen. In the lateral direction, the LVDT cores 
were spring forced against teflon discs fixed to steel pins 
grouted at an average depth of 35 mm into the specimen. 
These discs were located midway along the length of the 
Sides at a height dependent upon the location of the radial 
extensemeters. 

In the longitudinal direction, displacements were 
measured by LVDTs located at each corner of the loading 
head. In addition, dial gauges were also located at each 


corner to provide quick visual monitoring during testing. 


3.2.2 Radial Strain 

The relative radial displacement was measured by LVDTs 
connected to extensometers located in small diameter drill 
holes adjacent to the tunnel. Bach extensometer consisted of 
a coaxial rod and tube grouted at their ends ina hole 


within the rock mass as shown in Figure 3.4. The relative 
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Figure 3.4 Typical Extensometer Set-up (to scale) 
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displacement of these two anchor points is measured by an 
LVDT fixed to the outer ends of ‘the coaxial rods; the inner 
rod is fixed to the core and the outer tube to the body of 
the LVDT. The average radial strain for each extensometer 
was calculated as the relative displacement divided by the 


gauge length (ie. the length between anchor center points). 


3.2.3 Tunnel Closure 

Tunnel closure measurements (i.e., the ratio of the 
change in tunnel diameter to the original diameter) were 
made along four diameters spaced 25 mm apart along the 
tunnel axis and 45° apart in the plane perpendicualr to the 
tunnel axis. MeaSurements were made by LVDTs mounted inside 
the tunnel on a stand secured to the sample base plate as 
Shown in Figure 3.5. Two LVDTs are mounted on the same axis, 
their cores pressed to the tunnel wall by a spring, to 


measure the closure of opposite points on the tunnel wall. 


3.2.4 Sample Instrumentation 

Deformation measuring instruments within the sample are 
located by reference to a cylindrical coordinate system. The 
vertical axis corresponds with the tunnel axis and the polar 
axis with the principal field stress direction. The top 
surface of the specimen represents Station 0.0 on the tunnel 
axis. Station numbers correspond to the vertical distance 
from the top surface in millimetres. The mid point of the 


gauge length is used in locating radial extensometers. 
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spring and guide 
assembly 


HP 24DCDT-100 
displacement transducer 


levelling screw 
securing screw 


magnet 


Figure 3.5 Tunnel Convergence Measurement Device (electrical 
wires omitted for clarity) 
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Plan and sectional views of the instrumentation layout 
employed in Sample MC-5 are presented in Figure 3.6. 

Lateral boundary displacements were measured at a 
radius of approximately 270 mm from the tunnel axis at 
Station 142.5 (L1 to L4/0néeFigure 336) . 

Radial strain measurements were made by a total of 16 
extensometers located in two rings concentric with the 
tunnel axis. The first ring (R1) was located at a radius of 
103.6+1.5 mm from the tunnel axis at Station 92.5 while the 
second ring (R2) was located at a radius of 140.9+2.3 mm at 
Station 117.5. Average gauge lengths were 32 and 47 mm 
respectively. 

For Tests MC-5.04 and MC-5.05 and from MC-5.10 on, 
tunnel closures parallel and perpendicular to the major 
joint set were measured at Station 120.6 and 44.4 
respectively. Closures parallel and perpendicular to the 
major principal stress direction were measured at 
Stations 95.2 and 69.8 respectively. For tests MC-5.06 to 
MC5.09 closures parallel and perpendicular to the major 
joint set were meaSured at Stations 44.4 and 120.6 
respectively while closures parallel and perpendicular to 
the major principal stress direction were measured at 


Stations 69.8 and 95.2 respectively. 


322,452 .caliplenMC=6 
Plan and sectional views of the instrumentation layout 


employed in Sample MC-6 are presented in Figure 3.7. 
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Instrumentation Layout - Sample MC-5 
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Figure 3.7 Internal Instrumentation Layout - Sample MC-6 
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Lateral boundary displacements were measured at a 

radius of approximately 275 mm from the tunnel axis at 
babtoneio3.0 CiiletoL4 "on Fagore 1377,)5 

Radial strain measurements were made by a total of 24 
extensometers located in three rings around the tunnel axis. 
The first ring (R1) of eight extensometers was located at a 
Gadius) of 9050£2.7 mm at Stammone7s.0, the second ring (R2) 
ateapladtus «Of 615. 3-228 mm at Station 103.0, and ithesthird 
ringp(R3) at a radius of 7i3949+34i\mm at Station 128.0. 
Average gauge lengths were 34, 33 and 33 mm respectively. 
Rings R1 and R2 were monitored through Test MC-6.04. From 
Test MC-6.05 on, rings R2 and R3 were monitored as ring R1 
had been excavated during tunnel widening. Two special 
extensometers (S1 and $2) were installed at Stations 55.0 
and 67.5 with gauge lengths corresponding to the initial 
tunnel diameter. These instruments were monitored until cut 
through during the drilling operation of Test MC-6.02. 

Closure measurements were made at Stations 63.5 and 
83.9 parallel and perpendicular to the major joint set and 
at Stations 114.3 and 38.1 parallel and perpendicular to the 


major principal stress idirection for all tests. 


3.2.40 Sample MC-/ 
Plan and sectional views of the instrumentation layout 
employed in Sample MC-7 are presented in Figure 3.8. 
Lateral boundary displacements were measured at a 
radius of approximately 257 mm from the Sample centre at 


Station 131 for Sides 1 and 4 (L1 and L4), at Station 181 
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Figure 3.8 Internal Instrumentation Layout - Sample MC-7 
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fannSidex2ath2) candvateStationhh3éebor £Side i3.¥(b3)h. 
Radial strain measurements were made by a total of 24 
extensometers arranged in the following manner: 
Station 56: r =©104.141.52 mm at-orientations of @=45°, 
COPri225% gandn270 & 
Stabionséticnm= 128u6+3<0tmmeat 45° intervalsniromed=0° 
to 3452 inclusive: 
Station 106: r = 104.2+1.6 mm at 45° intervals from 6=0° 


£OgShS°? inclusive; 


Station 131: r 128.0 mm at an orientation of 6=90°: 
and 


Station 1563: r 


103.741.9 mm at orientations of 

G2450 5,902 ¥anda225 ci 
where r is the distance from the centre of the sample to the 
mid point of the gauge length. The average gauge length was 
35umm. 

Tunnel closures parallel and perpendicular to the major 
joint set were measured at Stations 113 and 88 and parallel 
and perpendicular to the major principal stress direction at 
Stations 139 and 62 respectively from Test MC-7.10 on. Prior 
to this the opening was not of sufficient diameter to allow 


closure measurements. 


3.3 Description of the Test Facility 
The equipment and methodology employed in this test 
program are essentially the same as used by Kaiser (1979), 


Guenot (1979) and Sargent (1981) with the exception of the 
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tunnel excavation procedure. The test facility was designed 

by Kaiser based on a similar facility developed by Heuer and 
Hendron (1971). A detailed description of the test facility 

is presented by Kaiser and Morgenstern (1981a). Only a brief 
description is repeated here. 

The test frame shown in Figure 3.9 was designed to 
permit long term testing of large specimens (up to 
approximately 600 x 600 x 200 mm) under plane strain 
conditions. In the lateral directions loads of up to 1300 kN 
are transferred from the hydraulic rams through load cells 
and load distribution triangles to the sample. Longitudinal 
expansion is nulled (thereby maintaining plane strain 
conditions) by applying load to the longitudinal loading 
head through 4 hydraulic rams. The hydraulic pressure system 
can maintain presSures over a period of several days with 
fluctuations of less than “one=percent. 

The advancement of a tunnel under stress was performed 
by full face drilling of the tunnel employing the equipment 
shown in Figure 3.10. Tunnels were advanced using diamond 
bit core barrels attached to a heavy duty Milwaukee 
DyNCdrTieOperating at /S0mrpm.sAlreatel008to. 150, kPaawas 
used to remove cuttings and to cool the bit. During the 
excavations performed on Sample MC-6 the drill was 
mechanically advanced. To provide better control over the 
rate of advance, the advance machanism was fitted with a 
Drayton RQ motor and appropriate gearing prior to drilling 


of Sample MC-7. Reaction for the drill was provided by rigid 
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Figure 3.9 Side View of Test Frame 
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Figure 3.10 Tunnel Excavation Equipment 
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connections between the drill and support shaft which was in 
turn hydraulically jacked against the concrete ceiling. 
Industrial vacuum cleaners were employed to collect the dust 
from the drilling operations and to clean the excavation. 
The data from all load cells and deformation recording 
devices are recorded automatically on magnetic tape and then 
transferred to the central computer. The data is then 
processed and the results plotted. The original test results 
are kept on file at the Department of Civil Engineering, 
Geotechnical Division, the University of Alberta by P.K. 


Kaiser. 


3.4 Loading Sequence 

The loading sequence of a given specimen consists ot 
two distinct phases. During the first phase, a series of 
tests is performed on the intact specimen to determine the 
deformation properties of the material and their spatial 
distribution. In the second phase the behaviour of the 
loaded specimens during the excavation or widening of a 
tunnel or the response of the specimen containing a tunnel 
to a change in load is monitored. 

The loading history of Samples MC-5, MC-6 and MC-7 are 
Mresented inser ouces 3. 1,1mto 3.13 respectavel ys In these 
figures the average vertical (i.e., major) field stress, the 
Pielaestress Latio(racio Of mnOrizontalmtosVenuicalst1e1d 
stress), the Poisson's ratio in the longitudinal direction 


(ratio of the stress in the longitudinal direction to the 
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Figure 3.11 Loading History of Sample MC-5 
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Figure 3.12 Loading History of Sample MC-6 
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Figure 3.13 Loading History of Sample MC-7 
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sum of the other two stresses) and the loading rate during 


each single loading stage are given as functions of time. 


4, DETERMINATION OF MATERIAL PROPERTIES 

Coal is the product of the physical and chemical 
alteration of vegetal matter accumulated in low energy 
non-marine to mixed marine environments. Because of the 
depositional environment coals commonly include layers or 
lenses of sandstones, siltstones or shales. The coal itself 
contains a variety of plant tissue (macerals) in different 
stages of preservation. This heterogeneity of the coal 
material will therefore result in non-uniform responses to 
stress change and thereby complicate the selection of 
representative properties. 

In the following sections the parameters were selected 
to describe the material property representing the average 
response; Significant local deviations from this mean 
response were observed. This particular investigation was 
concerned with the deformation properties. The strength 
properties reported summarize the work of earlier 


investigators. 


4.1 Strength Properties 

Investigations of the strength of the coal were 
conducted by Noonan (1972), Kaiser (1979) and Guenot (1979). 
Both direct shear and triaxial compression tests were 
employed. The unconfined compression strength of the coal 
was found to lie between 8 and 12 MPa. A bi-linear Coulomb 
failure criterion best described the strength under a 


multi-axial stress state. At normal stresses below 1 MPa, 
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the strength parameters, c, the cohesion intercept, and #¢, 
the angle of internal friction were approximately 0 MPa and 
60-70° respectively. For normal stresses above 1 MPa, a 
cohesion intercept of 2-3 MPa and an angle of internal 
friction of 30° were found to best fit the experimental 
results. The cohesion intercept varied inversely with the 
contamuity, ofmmointinget( iLie.@ipthesmatiosofsithesareal of the 
open joint surface to the total area of the plane containing 
the joint). The ultimate friction angle was found to be 30°. 
The axial strain at the yield point was 0.1-0.2 % at 
confining stresses less than 1 MPa rising to 0.5 % at higher 
confining pressures. Axial failure strains varied from 0.5 
to 1.6 % at low confining stresses. Under higher confinement 
failure strains up to 3.75 % were recorded. The strain 
values presented above were obtained from triaxial 
compression tests conducted at constant strain rates of 0.3 


toi Os 8 X%/nri. 


4,2 Deformation Properties 

In order to ascertain the constitutive relationship 
between applied stress and resultant strain for this coal, a 
series of tests were performed on each sample prior to 
tunnel excavation. Parameters describing both the 
time-independent behavior and time-dependent deformation 


behavior were thus determined. 
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4.2.1 Time-Independent Behavior 

The observed instantaneous or time-independent response 
of the intact sample to an applied increment of stress was 
typical of a fissured rock. The stress-strain curves showed 
an initial upward concavity at low stresses (below 3 MPa on 
average) becoming linear at higher stresses (maximum applied 
stress: 12.5 MPa). First loading produced appreciable 
Strains as a result of joint and crack closure. Subsequent 
loadings produced significanly less straining and more 
repeatable stress-strain responses. These phenomena, as 
observed in the testing of Sample MC-7 are demonstrated in 
Figure 4.1. 

Triaxial compression tests at confining pressures up to 
10 MPa conducted by Kaiser (1979) and Guenot (1979) have 
shown that the coal is an elastic, brittle plastic material 
with a strain weakening eine Pisce Non-linearity, 
aside from the initial, is only observed at stresses 
approaching the peak strength of the coal or if the strain 
rate is changed. The linear portion of the stress-strain 
curve was employed in the determination of the elastic 
parameters. 

In order to determine representative parameters, a 
procedure such as that outlined by Ko and Gerstle (1976) in 
which all the coefficients of the compliance matrix are 
determined, should be followed. Howvever, such an 
investigation was beyond the scope of this research. In its 


place, the form of the constitutive elastic relationship was 
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Figure 4.1 Average Extensometer Response During Hydrostatic 
(N=1) Loading of the Intact Specimen - Sample MC-7 
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inferred from the sample response and appropriate parameters 
determined. 

The spatial arrangement of the radial extensometers 
allowed for the determination of any directional dependence 
of the deformation properties in the plane of bedding. Under 
conditions of equal lateral stress, no significant variation 
in extensometer response was observed in Samples MC-5 and 
MC-7 as shown in Table 4.1. Sample MC-6, however, 
demonstrated a marked orientation dependence being much 
stiffer in the 0° - 180° orientation. This higher stiffness 
was caused, most probably, by the clay-shale inclusion 
previously noted and as such did not reflect an inherent 
anisotropy in the bedding plane of the coal. Hence, for the 
purposes of this research, it was decided to consider the 
Sample to be isotropic is the plane of bedding. This allowed 
Simplification of the analysis of the time-independent 
response to isotropic or transverse-isotropic elastic 


formulations. 


4.2.2 Material Properties for Isotropic Elastic Formulation 
Describing the linear elasic, isotropic behavior of the 
coal required the determination of two parameters: Young's 
modulus, E, and Poisson's ratio, UV. Formulae for their 
determination are presented in Appendix A. 
Figure 4.2 shows the spatial distribution of Young's 
modulus, normalized to the mean value, as determined by 


Method 1 of Appendix A at a field stress of 7 MPa for 
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Table 4.1 Orientation Dependence of Extensometer Response 


SLOPE RATIO(1) 
ORIENTATION | Test MC-5.03(2) | Test MC-6.02(3) | Test MC-7.08(3) 


0 - 180 
45 - 225 
90 - 270 
JEKY G EAB) 


(1) Average slope of stress-strain curve for that orientation 
divided by the average slope of all orientations 


(2) Slopes determined from linear regression of stress-strain 
coordinates at the beginning of each creep stage in 
multi-stage creep tests 


(3) Slopes determined from tangent to stress-strain curves 
at a field stress of 7 MPa 
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Figure 4.2 Spatial Distribution of Young's Modulus : 
a) Sample MC-5 b) Sample MC-6 and c) Sample MC-7 
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Samples MC-5 to MC-7. Mean values for Young's modulus and 
Poisson's ratio for these specimens were 1.73 GPa, 0.34, 
1.40 GPa, 0.29 and 1.50 GPa, 0.23 respectively. Application 
of Method 2 to Sample MC-7 (only intact specimen tested at 
different N values) yielded values of Young's modulus and 
Poisson's ratio of 1.46 GPa and 0.23 which were very similar 
to those obtained by Method 1. The spatial variation in 
deformation properties can be related to material 
heterogeneities as mentioned previously. 

The shape of the stress-strain curve and hence the 
derived elastic parameters depends, under isothermal 
conditions, on the combined effects of stress level and 
loading rate. In an attempt to ascertain the relative 
influence of the stress level, the elastic parameters were 
determined at stress increments of 1 MPa from an initial 
value of 3 MPa. 

Figures 4.3 to 4.5 show the mean and range of values 
for the elastic parameters determined by Method 1 for 
Samples MC-5 to MC-7 respectively. All specimens exhibit an 
increase in Poisson's ratio with increased stress. This 
phenomenon may be associated with widening of the joint 
space [Lama and Vutukuri (1978)] or dilatancy [Goodman and 
Dubois (1972)]. The opposite effect was noted when Method 2 
was applied for Sample MC-7 as shown in Figure 4.6. The 
effects of joint spacing.,and dilatancy are more suppressed 
Gugingutesting Of an intact ispecimen. Increased i stress may 


result in compression of the pores and increased crack 
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Figure 4,3 Variation in Elastic Parameters with Stress - 
Sample MC-5, Method 1 
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Figure 4.4 Variation in Elastic Parameters with Stress - 
Sample MC-6, Method 1 
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Figure 4.5 Variation in Elastic Parameters with Stress - 
Sample MC-7, Method 1 
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Figure 4.6 Variation in Elastic Parameters with Stress - 
Sample MC-7, Method 2 
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closuresthereby resulting in the lower- values: for Poisson's 
ratio noted. 

The influence of stress level on Young's modulus was 
inconsistent. Sample MC-5 (Figure 4.3) exhibited a decrease 
in modulus with increasing stress while Samples MC-6 and 
MC-7 (Figures 4.4 and 4.5) exhibited little change 
(employing Method 1). Application of Method 2 to the results 
from Sample MC-7 showed a direct correlation between modulus 
and stress as would normally be expected (Figure 4.6). 

While the results of this investigation regarding the 
influence of stress level on the elastic parameters are 
somewhat inconclusive it 1S most probable that the effects 
observed in Figure 4.6 are representative. This is due to 
the uncertainty of the assumed elastic stress distribution 
around the opening in Figures 4.3 to 4.5. Yielding of the 
near tunnel material, a distinct possibility, will certainly 
also influence the results presented in these figures. 

Results from investigations of the time-dependent 
response of the coal have shown that the magnitude of the 
creep rate increases rapidly at high stress levels (see 
Section 4.2.2). Mechanical constraints as well as test 
requirements restrict stress appliction to a maximum of 
approximately 0.005 MPa/sec. This is significantly below the 
rate of 0.07 - 0.7 MPA/sec commonly applied in static 
testing [Obert and Duvall (1967)]. Because of this slow 
loading rate, the strains recorded at higher stress levels 


will have increasing creep components which can result ina 
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non-linear, concave downwards stress-strain curve. This in 

turn results in reduced values of Young's modulus and 

possibly an increase in Poisson's ratio as the ratio of 
lateral to longitudinal creep strains may be greater than 
the ratio of elastic strains. Quantification of this effect 
was not possible given the test constraints, however, it is 
not expected to be any more significant that that due to the 

Stress level for the given rate of stress application. 

The effects of loading rate and stress level upon the 
elastic properties of the coal are to a certain extent, 
counterbalancing. It is therefore probable that the 
Parametensppresentedsane,neépresentatives of the material: 
This assumption is verified by two observations: 

1. The sample to sample variation in material properties is 
of the same order as those effects due to stress level 
and loading rate. 

2. The parameters determined for this coal are similar to 
those determined by others for various coals using 
Standard methods (see Table 4.2). 

Based on the above results, it may be concluded that 
under the assumption of linear elastic, isotropic behaviour, 
the coal has a Young's modulus of approximately 1.6 GPa and 
a Poisson's ratio of 0.3. However, in the longitudinal 
direction the stress required to maintain plane strain 
conditions is much less than predicted based on these 
parameters. For this reason, a linear elastic, Cransverse 


isotropic formulation is more probable. 
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Table 4.2 Elastic Properties of Coal 


Transverse 
Transverse 
Transverse 
Isotropic 
Isotropic 
Isotropic 
Isotropic 
Isotropic 
Transverse 
Transverse 
Transverse 
Transverse 
Transverse 
Isotropic 
Isotropic 
Isotropic 
Isotropic 
Isotropic 
Isotropic 
Isotropic 
Isotropic 
Or thotropic 
Or thotropic 
Or thotropic 
Or thotropic 
Or thotropic 


COAL | ELASTIC OUNG’ S MODULUS POISSON'S RATIO SHEAR MODULUS 
FORMULATION | Ex |Ey | Ez | xy | yx | xz [| zx | yz | zy | yz | xz / xy 
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NOTE: ALL MODULI IN GPa 


1 Dynamic, plane parallel plates 5x40 mm in X-section, 
Terms 3 to 8, 10 and 11 Coalification Series of Vitrain 
Schuyer et al (1954) 
Dynamic, 3.2x3.2x20-30 mm, Barnsley Hards, Morgans & Terry (1958) 
Static, 12.7 mm cube, Barnsley Hards, Morgans & Terry (1958) 
Static, 12.7 mm cube, Anthracite, Morgans & Terry (1958) 
Static, 38.1 mm cube, Barnsley Hards, Morgans & Terry (1958) 
Static, 38.1 mm cube, Anthracite, Morgans & terry (1958) 
Static, 25 mm dia. core, Bieniawski (1967) 
Static, 41 mm dia. core, Bieniawski (1967) 
Static, 50 mm cube, Bieniawski (1967) 
Static, 600x600x600-1200 mm, Bieniawski (1967) 
Static, 900x900x600-900 mm, Bieniawski (1967) 
Static, 1200x1200x600-1200 mm, Bieniawski (1967) 
Static, 1500x1500x600-1500 mm, Bieniawski (1967) 
14 Static, 100 mm cube, Pittsburg, Hawk & Ko (1980) 
15 Static, 100 mm , Old Ben, Hawk & Ko (1980) 
16 Static, 100 mm cube, York Canyon, Hawk & Ko (1980) 
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17 Static, 100 Beehive, Hawk & Ko (1980) 
18 Static, 100 Paonia, Hawk & Ko (1980) 
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4.2.3 Material Properties for Transverse Isotropic 
Formulation 

A total of five parameters are required to fully 
describe the linear elastic, transverse isotropic behaviour 
of a material: 

1. Young's modulus in the plane of isotropy (E ); 
2. Young's modulus normal to the plane of Meera (E); 
3. Poisson's ratio relating strains in the plane of ; 
isotropy (VU ); 
ré 
4. Poisson's ratio relating the strain in the plane of 
WEOELOpyecOu the Strainenormal tomity WU: 
5. Shear modulus relating shear stresses fi shear strains 
out of the plane of isotropy (G ). 
The full constitutive relationship aap the method of 
determining the five constants are presented in Appendix A. 
Data sufficient for the determination of these parameters 
was only available from the results of the tests performed 
on Sample MC-7. 

A hydrostatic compression test (all three principal 
Stresses equal; Test MC-7.03) enabled determination of E as 
0.314 GPa. This most probably consititutes a lower hones) se 
this modulus because it was essentially an initial loading 
modulus. The reloading modulus could be significanly higher 
as previously demonstrated (see Figure 4.1). Poisson's ratio 
relating the in-plane strain to the normal strain, VU , was 

zr 


found to be 0.082+0.015. The above two values were then 


employed to calculated the remaining constants presented in 
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Table 4.3. Insufficient data was available to enable 


calculation of G-. 
Z 
Finite element simulation of a tunnel in a medium with 


the above properties (assumed shear modulus equivalent to 
that in an isotropic medium) has shown that under plane 
Strain conditions, the stresses and displacements around the 
opening are not significantly different from those in an 
isotropic medium with the properties as given in the 
previous section. For this reason, subsequent analyses of 


test results assume linear elastic, isotropic behavior. 


4.2.4 Time-Dependent Behavior 

The time-dependent deformation of a material under a 
constant stress (creep) iS a major engineering concern, 
particularly in underground works. In order to address this 
concern it 1S necessary to determine ‘a priori' the 
time-Cependent behaviour of the material. This is very 
difficult as the long term response of the material may not 
be predictable from the short term laboratory response. 

An idealized strain versus time curve for a material 
with time dependent strength and deformation properties at a 
constant stress 1s shown in Figure 4.7. This curve can ibe 
divided into four regions: 
1. Blastic response of medium to change in stress; 
2. Primary creep stage characterized by a continual 

decrease in the rate of straining; 


3. Secondary creep stage characterized by a constant strain 
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Table 4.3 Parameters from Linear Elastic, Transverse 


Tests Employed 


Yzr 
Ez 


BOGE Gm 20D 
.04 & 7.05 
FOB EG. 7.05 
SO6NG. 1,07 


= 0.082 
= 0.314 GPa 


Isotropic Formulation 


Poisson’s Ratio 


Young’ s Modulus 


er 
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STRAIN 


REGION: 


1. Elastic Strain 
2. Primary Creep 
3. Secondary Creep 
4. Tertiary Creep 


TIME 


Figure 4.7 Development of Strain with Time 
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rate. This secondary creep stage may be very short in 
durationsor may not exist at aldeinea brittle rock like 
coal; and 
4, Tertiary creep stage characterized by an increasing 
Strain rate resulting in the ultimate failure of the 
Specimen. 
The same three creep stages (Regions 2 - 4) as characterized 
above are also shown in a double logarithmic diagram of 
strain rate versus time in Figure 4.8. The primary and 
tertiary creep stages are shown as linear, however, this 
need not be the case. 

In order to define the creep response of the material 
in a simple manner, an empirical approach, based upon curve 
fitting of experimental data, is generally employed. In most 
rock engineering problems the applied loads are sufficiently 
low that the tertiary stage doesn't develop and hence, is 
neglected [Schwartz and Kolluru (1981)]. The existence of a 
secondary creep stage has not been rigorously proven [da 
Fontoura (1980)], although recent studies by Tocher and Ko 
(1980) on coal have shown that a region of constant strain 
rate is achieved after approximately 2000 hours. No evidence 
of a secondary creep stage was found for the coal of this 
study under testing of up to 120 hours duration [da Fontoura 
and Morgenstern (1981)]. As no testing performed during this 
study was of longer duration expressions for the primary 


creep of this coal will be sufficient. 
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Figure 4.8 Idealized Double Logarithmic Plot of Strain Rate 
vs Time 
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The time laws describing primary creep are generally of 
the exponential or power law form. Cruden (1971a) 
demonstrated the superiority of the power law for primary 
creep in rock and da Fontoura (1980) showed its suitability 
for this coal. The parameters required for this formulation 
are eaSily obtained from a log strain rate versus log time 
plot as shown in Figure 4.9. Expressions obtained from 
Strain rate vs. time data are inherently better than those 
obtained from Strain vs. time data [Cruden (1970)] and 
therefore this procedure was employed in the determination 
of the appropriate parameters. 

Most studies into the creep of materials under constant 
Stress have been concerned with the creep resulting from a 
deviatoric stress condition (deviatoric creep). However, in 
this study, in order to describe the time dependent 
deformations observed in the laboratory it was necessary to 
develope expressions for the creep resulting from a decrease 
in stress (stress relief creep or relaxation) and that 
resulting from the uniform compressive stress (hydrostatic 
creep). For the purposes of this research it was assumed 


that direct superposition of all components was admissable. 


4.2.4.1 Hydrostatic Creep 

The term ‘hydrostatic creep’ in this study refers to 
those time-dependent deformations resulting from an equal 
Paterapestress, NOwLOnG tudinal straimecond ree oneandenctsa 
true hydrostatic stress condition..A linear relationship 


between the log of the strain rate and the log of time was 
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Figure 4.9 Determination of Respresentative Time-Dependent 
Parameters 


a2 


observed for each extensometer during single increment or 
the initial stage of multi-stage creep tests performed on 
the intact specimen under these conditions. Because of this 
a power law form of the 'hydrostatic' creep expression was 
adopted. 

Each of the extensometer responses was analyzed and a 
set of 'a' and 'm' parameters, representing the mean of the 


values, selected for a particular test. The hardening 


5 ‘ 


parameter 'm' which represents the rate of dissipation of 
strain rate with time was found to be relatively invariant 
of the 'hydrostatic'’ stress level as shown in Figure 4.10. A 
value of 'm'=0.83 would seem appropriate for this material. 
The magnitude of the strain rate at unit time (one hour 
in this study) is defined as the parameter ‘a' and it 
represents the strain potential. This parameter exhibited a 
strong stress dependency as shown in pices 4.11. Least 
square regression analysis of the data resulted in the three 
expressions for the 'a’ parameter shown. All had very 
Similar correlation coefficients (r?=0.68) and would yield 
Similar results in the range of interest (3 to 12.5 MPa). 
However, it was decided to adopt the power law relationship 
for the following reasons: 
1. The exponential relationship would tend to overestimate 
the parameter in both the low and high stress ranges. 
Indeed its use would result in predictions of 


hydrostatic creep at zero stress which is a physical 


impossibility under conditions of constant temperature 
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Figure 4.10 Variation in 'm' Parameter with Increasing 
Stress 


73 


i 


, ne 
i. a ciel = “i 


aa wt Teese, 5 


oo 


I 
! 
| 
I 
[ 
! 
3 
~ 
Cc 
co) 
Cc 
(o) 
oO 
«< 
lJ 


Linear 


=) 
(oe) wm co) 3 (o) 
N = = Oo 
° = S im S 
o oO ro) ro) 2) 


= 


AY/%,) YSLAWVYVd .D, 


Figure 4.11 Variation in 'a' Parameter with Increasing 
Stress 
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and humidity. 

2. The linear relationship suggested the existence of a 
threshold stress level beneath which no hydrostatic 
creep would occur. The existance of such a level is 
subject to question. 

Based on the above results the expression selected to 
describe the hydrostatic creep rate component was: 
ig ag 


eo td) BAG be Eqn. 4.1 
h h 
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Ao = hydrostatic stress increment in MPa 
12 : = 1.439 

‘e = time in hours 

n = 0.83. 


4.2.4.2 Stress Relief Creep 

The recovery of the sample upon removal of load was 
utilized to determine an appropriate expression for the 
stress relief creep. Again a linear relationship was 


observed in the double logarithmic plots of strain rate 


' ' 


versus time allowing determination of representative ‘a' and 
‘m' parameters. 

The 'm' parameter wasS again relatively constant at a 
Value of 1.05 as showneim Figure: 4.12. Schwartz and Kolluru 
(1981) reasoned that values for 'm' greater than 1 are 
physical impossibilities as the integration of such an 


expression would result in a creep function in which total 
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Figure 4.12 Variation in 'm' Parameter with Decreasing 
Stress 
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Strains decrease with time. However, this only is true if 
ene considers the limits of integration as being from t)=0 
SO some stime t2. Cruden (1971b\asuqgests that the east 
value of t; is the time taken for the stress wave to travel 
through the specimen (approximately equal to 2 x 10-* hours 
if the stress wave travels at about the speed of sound). Use 
of this lower limit for t, results in an expression in which 
creep strains increase with time. 

The 'a' parameter exhibited a dependence on the 
magnitude of the stress relief as shown in Figure 4.13. 
Least square regression analysis of the data resulted in 
three expressions with similar but poor correlation 
eoefficients (r7*=0.14): Cruden (1971by in a study of the 
recovery of Pennant Sandstone from uniaxial compressive load 
found the recovery rates to be stress independent. He 
Surmised that normal recovery results fron motion along 
cracks where the friction has been reduced with time due to 
corrosion aided failure of asperities. As the load on the 
specimen during recovery was much less than the applied 
loads the recovery rate would only depend on the cracks 
closed at this low load and hence be independent of the 
maximum load (cracks closed at higher loads would only open 
upon load reduction). While Cruden's results and rationale 
are compelling, the data presented in Figure 4.13 do 
indicate a slight stress dependency. Hence, it was decided 
to adopt the power law relationship (more reasonable at low 


magnitudes of stress relief). 
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decrease in stress is: 
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Ao = stress decrease in Mpa; 
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S = 0.41: and 

p =017205 


4.2.4.3 Deviatoric Creep 

The creep of the coal under deviatoric loading was 
investigated by da Fontoura (1980). He concluded that the 
hardening parameter 'm' was on average equal to 0.9 and 
independent of the stress level. The stress level was 
defined as the ratio of the current stress to the current 
Strength, both referring to the same state of confinement 
(to determine the current strength a cohesion intercept of 
2 MPa and an angle of internal friction of 50° were selected 
as representing average conditions). He also concluded that 
the 'a' parameter could best be represented as an 
exponential function of the stress level. The resultant 
expression for the deviatoric creep rate was: 


Bo -q 
e*(t) = Ae t Eqn. 4.3 
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o stress level; and 
gq = 0.9. 

Only one single stage creep test of the intact specimen 
under deviatoric loading conditions (Test MC-7.05) was 
undertaken. The mean and range of the measured values of the 
radial strain rates both parallel and perpendicular to the 
major principal stress direction for this test are presented 
in Figure 4.14. 

Two methods, each involving different assumptions were 
employed in the derivation of an expression for the 
deviatoric creep rate of Sample MC-7 at the stress level 
(Ows21) fof Test MC-7.05. 

In the first method it was assumed that Equation 4.1 
correctly described the hydrostatic creep rate component. 
The deviatoric component could therefore be calculated as 
the difference between the mean measured creep rate parallel 
to the major principal stress direction and the predicted 
hydrostatic creep rate. The deviatoric creep rate thus 
calculated is shown in Figure 4.15. Least squares linear 
regression analysis yielded the following expression for 
Lhis curve: 
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which is approximately 3.2 times greater than the expression 


proposed by da Fontoura. 
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Figure 4.14 Mean and Range of Measured Creep Rates in 
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In the second method it was assumed that there was no 
volume change due to deviatoric creep. Therefore the 
deviatoric creep rate components parallel and perpendicular 
to the major principal stress direction would be equal but 
of opposite sign. The deviatoric component of the creep rate 
can then be calculated as one half the difference between 
the mean creep rates measured parallel and perpendicular to 
the major principal stress direction. No assumptions 
regarding the hydrostatic component are required. The 
resultant values of the deviatoric creep rate components are 
also shown in Figure 4.15. Least squares linear regression 
analysis yielded the following expression for this curve: 


=0 91 
el(t)s= 4268@2tx5 40° *t (47nr) Eqn 4.5 
d 


which is approximately 1.6 times BEE than the expression 
proposed by da Fontoura. 

No other testing of the intact Specimen at different 
stress levels was undertaken. Therefore a general expression 
for a stress level dependent deviatoric creep rate could not 
be established for Sample MC-7. 

No further refinement of the deviatoric creep rate 
expression was attempted and Equation 4.3 with the listed 


Parameters was assumed applicable. 


4.2.4.4 Multi-Stage Testing 
The expressions describing the creep rate components 


developed in the preceding sections are not directly 
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applicable to multi-stage testing (i.e., tests where the 


load is applied incrementally with large time intervals 


between successive increments). In multi-stage tests the 


representative components are assumed to superpose in the 


following manner: 


fe =—n =i 
e' (tp =eBAcoritet to) BAG et Eqn 4.6 
h h hy 
Ss 6) S =p 
et (t) =aGAcait izatt) +eGAcprt Eqn 4.7 
sr d oe 
OsF =e} ByaNes ke 
enact) = Ae ¢t + Ae t Eqn 4.8 
d 
where Ao = increment of hydrostatic stress applied at 
hy 
time to; 
Ao = decrement in stress achieved at time to; 
dq; 
and 
Ao = increment in stress level applied at time 
Cow 


Similar expressions are assumed applicable for successive 


stages. 


Recent studies by Cruden (1983) suggest that a more 


appropriate form of the incremental creep rate is given by: 
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So = initial uniaxial compressive stress; 
S; = higher uniaxial compressive stress; 
to = duration of creep under So; 
t = time of creep under S,; and 


n,m = material constants which determine the 
exponents of both the time and stress dependence 
of creep rate. 

An investigation of the suitability of this expression for 
the results recorded during this test program was not 


undertaken, however it does warrant further investigation. 


4.2.5 Verification of Time Dependent Relationships 

As the expressions describing the three creep rate 
components were independently determined it was decided to 
verify their general applicability by comparing their 
predicted response with that measured during both single and 
multi-stage external loading tests under hydrostatic and 
non-hydrostatic loading conditions. The predictions of creep 
rate were made assuming stresses and stress changes 
corresponding to stress distributions ina linear elastic, 
isotropic medium with the properties previously determined. 
No verification of the stress relief creep rate was 
attempted. 

The expression for the time-dependent response of the 
rock mass to the external (i.e., boundary) loading of the 
specimen containing a tunnel would be a combination of 


hydrostatic and deviatoric creep. The hydrostatic creep 
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component is compressional whereas the deviatoric component 
in’the radia] direction is extensional.» 1fisuperposition is 
acceptable the form of the radial creep rate expression 


would be: 


eli (ti) Ge teemtty) amar eats) Eqn. 4.9 


which when expanded yields: 


Lee) DO aC 
cut) e-eBAg at ——Sem=t Eqn. 410 
is h 

Comparison of the radial Strain rates predicted by 
Equation 4.10 using the previously determined parameters 
with those measured during external loading of the three 
specimens containing a 152 mm diameter tunnel are presented 
in Figures 4.1i6a to c. The shaded area in each figure 
represents the range of strain rates recorded (all 
extensional) by the ring of extensometers closest to the 
tunnel wall. 

The measured radial strain rates are well predicted for 
Sample MC-5, however, they are significantly underestimated 
for Samples MC-6 and MC-7 (in Sample MC-6 even the sense of 
the strain rate is incorrectly predicted). In these two 
specimens the error is believed to lie predominantly in the 
deviatoric creep rate component. Earlier it had been 
demonstrated that the deviatoric strain rate determined by 


da Fontoura's relationship (Equation 4.3) significantly 
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underestimated the rate observed in Sample MC-7 (by a factor 
of 3.2). Assuming the ratio of the actual deviatoric strain 
rate to that predicted by Equation 4.3 to be stress level 
independent one can arrive at a new prediction for the net 
radial strain rate for Sample MC-7 as shown by the dashed 
line in Figure 4.16c. This clearly demonstrates that the 
error lay in the deviatoric component. These results imply 
that the actual deviatoric rate in Sample MC-6 must be 
approximately five times that predicted by Equation 4.3. 

Figures 4.17 and 4.18 present a comparison of predicted 
(solid line) and observed (shaded area) radial strain rates 
under multi-stage testing of the intact specimens under 
‘hydrostatic! stress conditions. Both figures show an 
initial slight overestimation of the rate with a continuing 
inclination to underestimation, particularly at times less 
than one hour, with successive daca tan wes Overall, however, 
the hydrostatic creep rate is suitably represented. 

Figures 4.19 and 4.20 present a comparison of the 
observed and predicted net creep rate under non-hydrostatic 
multi-stage testing conditions. Also shown on each figure is 
the predicted hydrostatic creep rate (dashed line). Parallel 
to the major principal stress direction (Figure 4.19a) there 
is reasonably good correlation between the predicted and 
observed. As had been the case under hydrostatic loading, 
there is an intial slight overestimation of the net rate 
followed by increasing underestimation particularly at times 


less than one hour. 
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Figure 4.17 Comparison of Observed and Predicted Radial 
Strain Rates Under Multi-Stage Hydrostatic Loading - 
Sample MC-5 
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Figure 4.18 Comparison of Observed and Predicted Radial 
Strain Rates Under Multi-Stage Hydrostatic Loading - 
Sample MC-7 
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Figure 4.19 Comparison of Observed and Predicted Radial 
Strain Rates Parallel to the Maximum Field Stress - 
Sample MC-7 
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Figure 4.20 Comparison of Observed and Predicted Radial 
Strain Rates Perpendicular to the Maximum Field Stress - 
Sample MC-7 
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Perpendicular to the major principal stress direction 
(Figure 4.20), there is virtually no correlation between 
predicted and observed strain rates. All observed rates 
(shaded area) were compressional whereas the predicted (not 
shown) were extensional. The observed rate exceeded even the 
predicted hydrostatic rate. The predicted net rates were 
determined using the parameters suggested by da Fontoura 
(1980) for the deviatoric rate component. These had 
previously been shown to significantly underestimate the 
deviatoric rate in this sample. Use of the sample specific 
parameters, however, would have resulted in even poorer 
predictions in this orientation. No explanation for this 
behaviour can be given at this time. 

The above results demonstrate that reasonably good 
approximations to material behavior can be achieved under 
Single stage external loading Pend aia by superposition of 
the hydrostatic and deviatoric creep rate components. Under 
multi-stage loading conditions the superposition of 
hydrostatic creep rate terms as given by Equation 4.6 yields 
sufficiently accurate predictions of sample response. The 
deviatoric component of creep rate is not well represented 
by Equation 4.7 in multi-stage testing. As no alternate 
relationship is currently available, Equation 4.7 will 
continue to be employed but with due consideration of its 
dOuUpeLULe Validity. 

Care must be employed in the general application of all 


the above relationships as sample to sample variations can 
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be considerable. There is an inherent risk involved in 
ascribing universal parameters to a highly heterogeneous 
material such as coal. However, it is important to point out 
that the proposed equations describe the time-dependent 
behaviour of the coal mass with all its joints and other 
discontinuities and hence, are of more practical value than 
results from small samples of 'intact' coal. 

The accurate representation of material behavior under 
complex loading conditions was not the purpose of this 
researehenit waskto adentifyethet parametersmsuffdacient ly 
accurately to provide a basis for the understanding of the 
time-dependent processes occurring near an advancing tunnel 
face. For this purpose the proposed expressions should be 


adequate. 
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9. RESULTS OF LABORATORY INVESTIGATION 


5.1 Introduction 

Current tunnelling technology is very much dependent on 
the analysis of deformations (closure, radial displacement) 
around the excavation. As explained in Chapter 2, these 
analyses are employed to establish the validity of the 
design and assess the stability of the excavation. While the 
magnitude of these deformations are important, particularly 
in reference to the serviceability of the opening, it is 
their evolution with time that is most beneficial to 
ascertain. 

The stand-up time and rate of tunnel cloSure are 
dominated by the excavation method and sequence, the type of 
Support and it's installation procedure. The time-dependent 
tunnel wall convergence is initially dominated by the 
time-dependent stress redistribution processes resulting 
from the advance of the tunnel face. Far behind the tunnel 
face the time-dependent behaviour is controlled largely by 
the creep properties of the rock mass. These time-dependent 
deformation mechanisms must be understood to achieve 
accurate evaluation of instrument response. 

In order to critically evaluate the current deformation 
monitoring practice, specifically the measure of tunnel 
convergence and the measure of radial displacement by 
borehole extensometers, a series of tests were performed on 


three laboratory specimens. Both external loading tests and 
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excavation Simulation tests were conducted. The results of 


this extensive laboratory investigation follow. 


5.2 Behaviour During External Loading Tests 

The term external loading test refers to all tests in 
which the load was applied to the specimen containing a 
pre-drilled tunnel. Loading was either quick (continuous 
increase from zero field stress to desired level) or slow 
(discontinuous increase from zero field stress with variable 
number of creep stages). Details of this loading sequence 
are presented in Chapter 3. The tunnel convergence and 
radial strain were monitored during both load application 
and each creep stage. These represent the time-independent 
and time-dependent responses respectively. The results of 
the external loading tests are presented according to these 


Giwrsions: 


5.2.1 Time-Independent Behaviour 

The material property distribution within a sample is 
very non-uniform, as demonstrated by the contour plots of 
the spatial distribution of Young's modulus presented in 
Chapter 4 (Figure 4.2). As a result of this spatial 
variation, the extensometer response recorded will be very 
much location dependent and wide variations in instrument 
response can be expected. Similar variations in the field 
must also be expected. The range of extensometer response 


under hydrostatic loading of the intact specimen is 
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presented in Figure 5.1 for Sample MC-7 (typical of the 
response recorded for all samples). In this and subsequent 
figures, the numbers in the lower right corner designate the 
location of the instrumentation along the tunnel axis. 
Compressive strains range from 0.6% (6 = 45°, Sta. 56) to 
1.9% (@ = 90°, Sta. 56) at a stress level of 10 MPa. This is 
much in excess of the predicted elastic strain of 0.46% 
obtained using the value of Young's modulus and Poisson's 
ratio suggested in Chapter 4. The principal reason for this 
excess deformation results from the extreme non-linearity of 
the stress-strain relationship at low stresses (due to crack 
closure under load application or opening under stress 
relaxation) which is not accounted for in the simple linear 
elastic relationships employed. From this preliminary 
investigation of the extensometer response, it is evident 
that it will be strongly dependent on Hora material 
properties, stress level and orientation with respect to 
intersecting discontinuities. The detection of abnormal 
behaviour from strain magnitudes recorded by extensometers 
willerbe: complicated by these influences. 

Figure 5.2 presents the tunnel convergence and radial 
strains recorded during hydrostatic loading of Sample MC-5 
containing a 120 mm diameter unlined tunnel (Test MC-5-04). 
Closure significantly in excess of the average is observed 
across the 0° - 180° diameter, although the radial strains 
do not indicate any abnormality in this orientation. They do 


however, indicate abnormal behaviour near the tunnel wall 
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Figure 5.1 Extensometer 
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Figure 5.2 Tunnel Closure and Extensometer Response :; 
Test MC-5.04 
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(Sta. 92.5) approximately normal to the maximum closure 
(Danieveularlyeat -6=-—45°) SAswicderevartaclonminmec losure 
response is noted during the next loading sequence (Test 
MC-5.05) as shown in Figure 5.3. Maximum closure is again 
across the 0° - 760° diameter but now ‘ae@distinct minimum 
closure is observed across the 135° - 315° diameter. The 
closure is recoverable in the latter direction, but 
permanent closures of 2.5% to 4% are observed in the other 
Orientations suggesting that yielding has occurred. These 
results are confirmed by the extensometer readings which 
indicate yielding near the tunnel wall (Sta. 92.5). Yielding 
processes appear concentrated in the 135° - 270°/315° 
Orientations. It is of particular interest to note that 
yielding is occurring in the direction of minimum tunnel 
convergence and almost normal to the maximum convergence. 
This contradicts the oft-held assumption of large dilation 
in the yield zone. 

In the above analysis of the observed responses, it was 
not so much the magnitude of the response that allowed 
interpretation as the difference when compared to previous 
measurements at that location or to other instrument 
behaviour. Similar comparative procedures must also be 
employed in practice. 

Figure 5.4 presents the tunnel convergence and radial 
Strains measured during hydrostatic quick loading of the 108 
mm diameter tunnel of Sample MC-7 (Test MC-7.11). Maximum 


closure was observed across the 135° - 315° diameter, 
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Figure 5.3 Tunnel Closure and Extensometer Response : 
Test MC-5.05 
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however, the variation in closure was not as severe as that 
observed in Sample MC-5 under similar conditions (Figure 
5.2). No abnormal extensometer responses were noted during 
loading, however, upon load removal, in the 135° orientation 
both near (Sta. 106) and far (Sta. 81) from the tunnel wall, 
permanent, extensional strains of up to 0.7% suggest that 
some Slip along jointing occurred. In this case the 
Orientation of maximum convergence corresponded to the zone 
Ofeactivity. 

The results of the slow, hydrostatic loading of the 152 
mm diameter tunnel of Sample MC-7 (Test MC-7.15) are 
presented in Figure 5.5. Closure is maximum in the 135° - 
315° orientation and minimum in the 45° - 225° orientation. 
Abnormal straining, indicative of yielding processes, is 
again predominant in the 135° orientation with some 
expansion of the yield zone to the 90° orientation noted 
near the tunnel wall (Sta. 56, 106 and 156). Again the 
maximum convergence corresponded to the active zone as 
denoted by the extensometers. 

Figure 5.6 presents the measurements recorded during 
non-hydrostatic (N=0.5) quick loading of the 152 mm diameter 
tunnel of Sample MC-7 (Test MC-7.16). The loading response 
of the convergence gauges corresponds to that expected ina 
homogeneous, elastic material with minimum closure normal to 
the principal stress direction (90° - 270° orientation) and 
Maximum parallel to it (0° - 180°) orientation. Upon 


unloading the closures in ‘all *but across the ‘0% = 180° 
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diameter are essentially recoverable. The measured radial 
Strains continued to indicate yielding in the 135° 
Oorlentation with expansion through 90° to 45° (Sta. 81). An 
interesting observation with regard to the strains recorded 
in the 45° orientation is that strains on the order of 2% 
(extension) were recorded far from the tunnel wall (Sta. 81) 
while much less straining was observed near the wall (Sta. 
56, 106 and 156). There is no reason to discard the reading 
at Sta. 81 as the instrument performed satisfactorily in 
earlier tests. Therefore, it must be concluded that both 
anchor points of the near tunnel extensometers are moving in 
unison as a result of deep seated deformation processes. 

As a result of these investigations, it can be further 
concluded that while tunnel convergence may give an 
indication of yielding or other stress redistribution 
processes, it does not provide any information regarding the 
type or location of the process. Extensometers can provide 
that information, however, they must be properly located and 


ofesurfrcient ‘extent. 


5.2.2 Time-Dependent Behaviour 

The absolute magnitude of tunnel closure or radial 
Strain cannot generally be determined in the field due to 
the impracticality of instrumentation ahead of the face in 
deep tunnels or the limited face access during advance. It 
has thus been found convenient to use the rate of 


deformation rather than the magnitude of deformation as an 
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indicator of tunnel performance (see Chapter 2). For this 
reason, interpretation of the results of this test program 
will be based on the time-dependent observations presented 
in double logarithmic plots of deformation rate versus time. 
As shown in the previous section and discussed by 
Kaiser (1981d) and Kaiser and Morgenstern (1982), tunnel 
convergence reflects the overall behaviour of the 
underground opening but must be combined with extensometer 
measurements to delineate the yield or failure mechanisms, 
their locale and extent. The above observations led to the 
selection of the closure rate as an indicator of tunnel 
behaviour with the extensometer response employed to define 


the behaviour. 


De2e2.) Critical Closure Rates 

The time-dependent response of the rock mass depends, 
to a large extent, on the intensity of the applied load. At 
low stress levels, terminating (primary) creep will lead to 
a finite creep strain whereas at high stress levels, a 
transition from primary to tertiary creep with increasing 
creep rates leading to failure will be observed. It would 
therefore seem likely that there is a strain rate boundary 
separating the terminating from the accelerating creep 
range. This critical creep rate function, CCR(t), is shown 
schematically in Figure 5.7a. 

This behaviour is typical of that encountered in 
uniaxial compression tests. Under load P,, the observed 


creep rate function falls below the critical creep rate 
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function and creep terminates. Under P2, the observed creep 
rate function is greater in magnitude than under P, but 
still less than the CCR(t) and the creep terminates. Under 
loads P; and P, the strength of the material is exceeded. 
The observed rates are greater in magnitude than the CCR(t) 
and accelerating creep resulting in failure occurs. 

In tunnelling situations, initiation of yielding will 
not necesSarily lead to collapse but to a change in the rate 
of deformation due to the change in the rock mass structure 
and hence, its deformation properties. The type of response 
expected in tunnelling situations is shown schematically in 
Figure 5.7b. Under a field stress S, only primary creep and 
no yielding occurs (below CCR(t)). Under field stress S, 
some yielding occurs as noted by the increase in deformation 
rate at point A. Stress is redistributed around the opening 
and a higher terminating creep rate curve is reached. This 
creep rate is called the post-yield creep rate, PCR(t). It 
is the rate that is reached after propogation has 
terminated. Besides other factors, it will depend on the 
extent and shape of the yield zone generated under field 
stress S2z. Under an even higher field stress S3 yielding may 
6cGaum fitustb@atipointeA andvagainvat point BeeAt this point 
the structure may collapse and become unstable if the degree 
of freedom is sufficiently high to create a collapse 
mechanism. Under S3; the creep rate function is such that no 
stable equilibrium is possible. The rate function at which 


thishunstablesiailure propagation occurs 1s called the 
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unstable propagation rate function, UPR(t). The UPR(t) can 
only be reached in near surface openings or in cohesionless 
materials where unstable propagation is possible. 

Figure 5.8 summarizes the observed closure rates 
observed in Test MC-4.2 [Guenot (1979)]. According to the 
rate function definitions given previously, two closure rate 
functions have been sketched. A linear CCR(t) was chosen to 
follow the upper limit of the non-yielding creep curves 
(shaded region). The PCR(t) was assumed to be parallel to 
the CCR(t) and was drawn touching the highest rates measured 
during the test at 14 MPa when stable yield zone propagation 
was observed. The UPR(t) could not be accurately located as 
no tests leading to collapse were executed. 

For simplicity and comparison with other test results, 
it was decided to consider rates at t = 1, 5 and 20 hours 
rather than the entire creep rate Function. The applicable 
values obtained from Figure 5.8 are: 

t=1 hr CCR(1)=30 ; PCR(1)=80-100 
t=5 hrs CGRG5) =9 Fee PGRC5) =25—30 XeLO ae ey iG 
t=20 hrs GCR(20)=3) 2" PER(20)=8-10 

The two creep rate functions defined in this manner 
allow characterization of the tunnel behaviour. 

The direct application of the closure rate functions, 
determined from external loading tests, in the evaluation of 
tunnel behaviour is questionable as the condition of 
external loading does not correspond with reality. However, 


if the processes are independentof the loading sequence, 
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Figure 5.8 Typical Tunnel Closure Rates During Test MC-4.2 
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then they are applicable. In order to determine the loading 
history dependence, the range of closure rate curves for 
tests with a stress ratio, N, of unity (shaded areas) were 
compared with the the mean compressional strain rates of the 
intact coal sample in Figure 5.9. 

Comparing pre- and post-excavation test results from 
Sample MC-5, it can be seen that the tunnel closure rates 
are consistently one order of magnitude (or more) greater 
than the mean compressional strain rate. Hence the component 
of pre-excavation creep in tests performed on Sample MC-5 is 
negligible and the measured rates correspond closely to the 
rates that would be measured in a test correctly modelling 
reality. 

Because of the higher creep potentials of Samples MC-6 
and particularly MC-7 under compression, a Significant 
portion of the tunnel closure rate eRe be due to the 
compression of the coal block in external loading tests. 
Comparison of closure rates measured during external loading 
and an excavation simulation test on Sample MC-7 (Figure 
5.10) however, indicates that there is little, if any 
difference and that the closure rate functions as defined 
may be applied. 

The evaluation of tunnel behaviour by closure rate 
measurements during excavation simulation and when the 
tunnel diameter was less than 108 mm was not possible. In 
these situations, tunnel behaviour was evaluated by radial 


creep strain rate functions defined in a manner analogous to 
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Figure 5.9 Compressional Strain Rates for Intact Samples and 
Ranges of Tunnel Closure Rates for Samples MC-5 to MC-7 
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Figure 5.10 Comparison of Observed Tunnel Closure Rates 
under Excavation (Test MC-7.10) and External Loading 
(Test MC-7.11) Conditions 


those presented for closure. In order to avoid confusion 


with the closure rate functions new labels have been 


assigned to the strain rate functions as follows: 


CSR(t 
PSR(t 


USR(t 


) corresponding to the CCR(t); 
) corresponding to the PCR(t); and 


) corresponding to the UPR(t). 


While the compressional creep component could be 


neglected in evaluating the tunnel closure rates, 


be disregarded in the evaluation of radial creep strain 


it cannot 


rates. The creep strain increment observed in the test may 


be dominated by the sample compression depending on the 


magnitudes of 


moduli[Kaiser 


the elastic (i.e., Young's) and bulk 


and Morgenstern (98:2) pleqnsai0): 


In an isotropic elastic medium under conditions of 


plane strain and equal lateral field stresses (N=1), the 


time-independent closure resulting from the excavation of a 


tunnel is given by: 


u(a)/a = Ppo/2G Eqn. 


where u(a)= 
a = 
Doma 
G 


and the time- 


€ 


tunnel wall displacement; 
tunnel radius; 

field stress; and 

shear modulus. 


independent radial strain by: 
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where r = distance from tunnel centre to the point of 
interest. 
Thesratro or, tunnel closure to radial stratn for thevelastic 


case 1s then given by: 


lula)/al/e = =(r/a)¢ EBqn.5.3 
3 
Assuming that the same relationship holds for the 
time-dependent behaviour under simulated field conditions, 
the appropriate values for the radial creep rate functions 
are -(r/a)? times the closure rate functions. For external 
loading tests the required values of the strain rate 


functions may be obtained from the following relationships: 


CSR i) Sie a /) 2 ex COR +) ata COMPiCt) 
PSRit) =) —(a/n)? x PCR(t). + COMP Gt) Eqn. 5.4 
USR(t) s=eatafr)? x UPR( tite COMP (Lt) 

where COMP(t) = the compressional strain rate at time (t) 


calculated assuming an elastic stress distribution. 


ene Cee varuallOny OL leStmReSUrts 

Tables 5.1 and 5.2 summarize the stress levels at which 
the CCR(20) or CSR(20) and PCR(20) or PSR(20) were reached 
for the mean (in brackets), the maximum and the minimum 
tunne hiclosure ‘or radial’ strainerate. These stressimlevels 
were determined from the plots of stress versus log closure 


rate presented in Appendix B. In the same tables are listed 
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Table 5.1 Predicted Stress Level for CCR(20) and PCR(20) and 
Observed Stress Level where Propogation Clearly Observed 
from Extensometers - Samples MC-5 and MC-6 


PREDICTED STRESS LEVEL 
CCR(20) PCR(20) 
) Mean ( ) Mean 


13.0-(>13.5)->13.5 


Observed Max. 
Radial Creep Strain 
ncrement (Extensometer) 


Predicted Observed 


<11.0-(11.5)->13.5 


<11.0->11.0 >11.0 


after unloading 


ee ee 
(AMBIes : ei ei poe] ee | ee 
pa pe ue 


15-<0.22 
13-<0.1% 


13-<0.1% (comp. ) ** 
13-<0.8% 


15-<0. 8% 
13-<0,.2% 
10-<0.1% 


15-<0.6% 
13-<0.2% 
10-<0.1% 


(MPa) {MPa} [loner ring strain [MPa - Mode] [MPa - Mode] 
recorded during entire 
Creep stage] 


[MPa - %] 


* only observed In last measurement; prior to this, the maximum creep strain Increment was 0. 13% Indicating 
Inttlation rather than propagation 

ae nearest strain measurement at 36 mm from tunnel wall; near tunnel! strain not observed, hence, ylelding 
May not be detected by Instruments. 
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the observed maximum radial creep strain increments 
occurring during the corresponding creep stage indicated by 
the stress level. The mode of behaviour was predicted by 
comparisons of the actual closure or radial strain rate with 
the proposed rate functions at a time of 20 hours. The 
observed mode of behaviour was determined from the creep 
strain increment in the following manner. If the creep 
Strain increment was less than 0.2%/20 hrs but obviously 
larger than at lower stress levels "initiation" (I) was 
assumed: if it was about 0.2 to 0.4%/20 hrs 
"initiation-propagation transition" (IP) was assigned; and 
1f ittmwasS greater than 0.6%/20 hrs "propagation yielding” 
(P) was assigned. This latter limit was selected 
arbitrarily, however, a creep strain of 0.6%/20 hrs will 
clearly lead to accumulated strains in excess of the failure 
Strain of coal and hence does eee wa paca: ac yielding. 
In the final column it is shown whether prediction and 
observation agree. 

Excellent agreement between prediction and observation 
were achieved for Samples MC-5 and MC-6. Good agreement was 
achieved for Sample MC-7 (Table 5.2) Apparent conflicts 
between predicted and observed behaviour for Sample MC-7 can 
be resolved by comparison of the full closure rate functions 
with the double log plots of closure rate versus time or 
consideration of test limitations as follows: 

Tests Mc-7.08 and MC-7.10 pre-widening - Because of future 


test requirements no instrumention was sufficiently close to 
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the tunnel wall to detect any yielding; 

Tests MC-7.13 pre-widening and MC-7.14 - While the limits of 

the appropriate functions are not exceeded at a time of 20 

hours, an examination of the complete curves shown in Figure 

5.11 show jumps indicative of yielding in the closure rate 

curves at earlier times. Two possible explanations for the 

failure of the proposed CCR in indicating the onset of 
yielding are: 

1. The values for the cnitical rate functions are too high 
for this specimen. Reduction of the functions by 
approximately 30% would yield more accurate predictions; 
and 

2. The behaviour mode is different than those previously 
described. Early high rates indicative of local slip or 
local yield are rapidly dissipated by stress 
redistribution without global yielding (i.e., the 
overall structure of the rock mass is not signifigantly 
altered). The longer term behaviour is that of a 
non-yielding tunnel even though the local strain 
magnitudes indicate that yielding has occurred. 

Test MC-7.15 - Examination of the complete curves show 

closure rates in excess of those defined by the CCR curve 

and approaching thegPCR curvejat a field stress of 12.5 MPa 

(Figure 5.12a). Figure 5.12b shows that at a stress of 14 

MPa the closure rate waS again in excess of the CCR curve 

and then exceeded the PCR curve at a time of 20 hours (this 


late surge in closure rate had been previously ascribed to 
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10-9 MC-7.13 
12.5MPa 


TUNNEL CLOSURE RATE (%/hr 


10' 10° aol 402 408 
TIME (hrs) 


MC-7.14 
12.5 MPa 


10' 10° of 402 © 493 
TIME (hrs) 


Figure 5.11 Observed Tunnel Closure Rates : a) Test MG=7rale 
pre-widening; b) Test MC-7.14 
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Figure 5.12 Observed Tunnel Closure Rates - Test MC-7.15 
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instrument fluctuation). Hence, agreement with the observed 
was achieved; and 

Test MC-7.17 - Examination of the complete curves (Figures 
5.13 a to c) shows that while the prescribed limits were not 
exceeded at a time of 20 hours, they had been at an earlier 
time (at principal field stresses of 10 and 12.5 MPa). At a 
maximum field stress of 7.5 MPa, the CCR curve was not 
exceeded, however a marked increase in closure rate between 
1 and 2 hours indicated that some yielding had occurred. 
Again, complete curves show agreement between predicted and 
observed behaviour even though at the specified time of 20 
hours, agreement had not been observed. 

It may therefore be concluded that the mode of 
behaviour could be predicted with reasonable accuracy for 
all external loading tests. The closure rate trend or change 
in cloSure rate seem to provide an excellent indication of 
the proximity and extent of stress redistribution processes 
(even though neither the failure mechanism nor its location 
can be determined) and hence is a good measure to evaluate 
the short and long-term stability of an underground opening. 

The derivation of displacement rate criteria for tunnel 
evaluation and their subsequent verification in model tests 
resulted in some important observations in the use of 
displacement rate criteria. These are as follows: 

1. The agreement achieved between predicted and observed 
behaviour in the three specimens indicate that the 


concept of critical rate functions should be suitable 
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Figure 5.13 Observed Tunnel Closure Rates - Test MGR ray 
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for stability evaluation practice; 
It is necessary to consider the complete (or at least a 
Significant portion of the) closure or strain rate curve 
and not just the rates at specific times. As shown in 
Figures 5.11 to 5.13 the observed rate may be less than 
the critical rates at the 'test' times but may have 
exceeded the limits at intermediate times. This in turn 
raises two further points: 
(i) Tunnel closure or radial strain measurement 
should be started as soon as possible after 
excavation and should be taken as often as possible 
to avoid masking of intermediate behaviour. In 
addition, reliable determination of tunnel advance 
effects are only possible if the early part of the 
rock response 1s known; and 
(ii) The tendency towards ividariteeeven (i.e. best 
fit) of the closure or strain rate versus time 
relationship in a double log plot should be 
discouraged. Apparent fluctuations in rate, that may 
be attributed to instrument error, may be indicative 
of real processes and have considerable bearing on 
the tunnel behaviour; 
Multi-directional tunnel closure or radial strain 
measurements are required to detect behaviour indicative 
of yielding. Too few readings may result in 
misinterpretation. Figures 5.11 to 5.13 show that often 


only one of the four diameters is experiencing closure 
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rates in excess of the prescribed limits. If measurement 
had not been made in this critical orientation the 
prediction of adverse behaviour could not have been 
made. This also supports the argument for local 
yielding; and 

4. Rates in excess of the prescribed closure or strain rate 
functions are not the exclusive indicator of behaviour 
mode. Anomalies in the double log plots of closure rate 
versus time as shown in Figure 5.13a at a time of 2 
hours are an indication of yielding or softening as 
well. It is the deviation from the normal trend that 
Sshouldebeyconsidered) rather than va (strict adherence to 
specific critical values. 

In practice, the problem of determining the site 
Specific closure or radial strain rate functions for a 
specific rock type remain to be ouye on However, by 
continuous monitoring of rates in stable, intermediate and 
unstable sections of an opening, one should be able to 
determine appropriate values. These can then be related to 
construction procedure and support capacity. This latter 
process requires the separation of the effects due to face 


advance from the time-dependent rock deformation. 


5.3 Behaviour During Excavation Simulation Tests 
A total of five tests on two specimens, Samples MC-6 
and MC-7, were performed to evaluate the behaviour of the 


tunnel during its advance. The processes involved could only 
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be determined from extensometer measurements as the 
recording of tunnel convergence during excavation was not 
possible. These extensometer measurements will to some 
extent be influenced during the excavation process by the 
proximity of the sample boundaries (perpendicular to the 
tunnel axis) that are one or less tunnel diameter from the 
meaSuring points. The effect of these boundaries in the 
direction of the tunnel axis should not dominate the 
behaviour because of the close simulation of the real 
condition by maintaining plane strain conditions in the 


plane perpendicular to the tunnel axis. 


5.3.1 Behaviour Near the Tunnel Face 

Panet and Guenot (1982) have analyzed the convergence 
behind the tunnel face and identified an apparently unique 
relationship between convergence related to tunnel advance 
and radius of the plastic zone. This relationship cannot be 
directly applied to the test results because of the lack of 
convergence meaSurements during excavation. Hence, it was 
necessary to determine a relationship between convergence 
and radial @strain. sThis#isitsimple?fon astunnelwintatlimear 
elastic material but complicated for tunnels in yielding 
LOCK. 

As a first approximation the test results were compared 
with the prediction for advance in a linear elastic 
material. The effect of time-dependent deformations are 


treated independently of the tunnel advance effect and are 
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assumed to superpose directly. 
For a circular tunnel in a linear elastic medium under 
uniform field stresses (N=1) the radial strains are given by 


Equation 5.2 and the displacements by: 
u(e)e =I po72C) Saez e) Ron. 5.5 


The ratio of radial strain to tunnel wall displacement, 


uta) cistthen® 
e (eu al G=R-a/r2 Delgts Shey 6 
r 


According to Panet and Guenot (1982) the tunnel wall 


displacement rate u'(x) for the elastic case is: 


ul (xo @=SE2u(@) (tka) 7A karx ye /ke ‘Helr 27 


where k OS 4m foundsby curves fitting) 

By = distance from tunnel face in the direction of 
the tunnel axis; and 

u(~)= ultimate elastic tunnel wall convergence 
(i.e., the difference between the convergence 
recorded when the face is at the measuring 
Station and when the face is at an infinite 
distance from the station). 


Employing the relationship defined by Equation 5.6 results 


in the following expression for the radial strain rate as a 
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Lunetnen cof ethe idrstance «Eromethestunnel face: 


e%(x)o=e2a2ulo)f (ka)/A( katx).]? /kr? Ean. °5.8 

is 
For the condition of tunnel widening, it is necessary to 
consider only the differential tunnel wall displacement, 
Au(@) (i.e., the difference in convergence, u(~), for the 
two tunnel sizes), due to incremental stress change to 
predict the rates. Equation 5.8 can be rewritten to yield an 
expression forjthedéradialistrainwrate asia vfunction of stime 


by use of the Chain Rule: 


ehGt jle=e(dae/dx )(dx/dt) EGnew .o 
r r 


iH 


where de /dx e'(x) as given by Equation 5.8; and 


ie r | 
adxu/dt senate ohetunnel vadvance, .d". 


Expanding Equation 5.9 and substituting d't for x yields: 


éheit)a= 2a? ulmddblk/(katd ste) J2ykn2 metry sah!) 
r 


A similar expression for the tunnel closure rate as a 


function of time can be derived as: 
u'(t) = -2u(@)d'[(ka)/(kat+td't)]*/ka eign, Seed 


The rates predicted by Equations 5.10 and 5.11 are 


applicable only during active advance of the tunnel. If 
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excavation or widening is interrupted or completed, the 

effect of advance is lost and only rates due to the 

time-dependent properties of the rock mass are measured. 

PnAordéreto 'verifymihe vahiditysofskquatiom 510 
comparisons of predicted response were made with the 
observed response in four excavation simulation tests, 
MG-~6502 -eMCo 208 MGS 270@dnaemMC-7.13.. Test MC-6.05 was not 
Similarly investigated due to the number of interruptions 
during excavation. Rather than just comparing the observed 
rate which included all time effects with the predicted rate 
due to face advance, it was decided to compare it with the 
predicted total response (obtained by incorporating those 
expressions developed in Chapter 4 for the time-dependent 
response of the rock mass). 

Rigorous Simulation of the actual response to 
excavation or widening was considered too refined for the 
purpose of this investigation, however, it was felt that 
reasonable approximations to the behaviour could be made by 
application of some simplifying assumptions: 

1. The material is homogeneous and isotropic with respect 
to both it's time-independent and time-dependent 
properties; 

2. The time-independent response is linearly elastic; 

3. The time-dependent response is as specified in Chapter 
4; 

4. The stress change due to the creation or widening of the 


tunnel at a given station along the tunnel axis occurs 
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immediately upon the arrival of the tunnel face. 
Comparison of this assumed stress response with observed 
stress changes during widening (Figure 5.14) as recorded 
by Korpach (1983) shows that this assumption is 
reasonable for preliminary analysis. Alternatively, a 
progressive change in stress such as proposed by Sakurai 

(1978) may be employed; and 
5. Direct superposition of the various components of 

time-dependent deformation is acceptable. 

Predictions of the radial strain rate were made, at a 
given station along the tunnel axis for the instruments 
nearest to and furthest from the tunnel wall (inner ring of 
extensometers only) in order to bound the expected 
behaviour. The predicted response was determined by 
Superposition of the following: 

1. Extensional strain rate due to #406 advance; 

2. Extensional strain rate due to stress relaxation; 

3. Extensional strain rate due to the deviatoric stress 
conditions created by the excavation and; 

4. Compressive stain rate due to the hydrostatic 
compression of the test specimen. 

The respective components along with the net rates after 

Superposition are shown in Figure 5.15 for Test MC-7.08 at 

Station 81, r/a=4.64 (extension negative). 

Figures 5.16 and 5.17 show the predicted (heavy solid 
line) and observed radial strain rates during and following 


excavation of a 120 mm diameter tunnel in Sample MC-6. While 
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Figure 5.14 Observed Stress Changes in Tunnel Widening 
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MC 7.08 
T/ OE 4 OGmola Gc | 


STRAIN RATE COMPONENTS 
——--—— Face Advance 
Deviatoric Stress 
——— Stress Relief 
—-— Hydrostatic Stress 


NET RATE 


RADIAL STRAIN RATE (°%/hr) 
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Figure 5.15 Predicted Radial Strain Rate Components - Tunnel 
Excavation (Test MC-7.08) 
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Figure 5.16 Comparison of Observed and Predicted Radial 
Strain Rates - Test MC-6.02, Parallel to Principal Stress 
Directions 
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Figure 5.17 Comparison to Observed and Predicted Radial 
Strain Rates - Test MC-6.02, 45° to Principal Stress 
Direction 


the post-excavation behavior is reasonably well predicted 
the strain rates during excavation are significantly 
overestimated. The poor agreement between predicted and 
observed radial strain rates is most probably a result of 
the boundary conditions of the model. Negligible shear 
stresses at the sample boundary result in little load 
transfer during advance to and slightly beyond the measuring 
Station. Strains and therefore strain rates are below those 
expected in an infinite medium. As the tunnel face 
approaches the other boundary which again has negligible 
shear stresses the load increases more rapidly than would be 
the case in reality and strains and strain rates exceed 
those predicted. These effects would be most pronounced 
during initial excavation. Behaviour such as this was 
observed during both excavation and widening operations. 
Figures 5.18 and 5.19 show the predicted (shaded area) 
and observed radial strain rates during and following 
excavation of the 42 mm diameter tunnel in Sample MC-7. The 
general trend of the radial strain rate is reasonably well 
predicted after a time of 0.2 hours. Prior to this the rates 
are overestimated as before. The face advance effect is 
clearly dominating the observed strain rates during 
excavation as demonstrated by the rapid increase in strain 
rate when drilling recommences following an interruption. 
Better correlation with the predicted value may have been 


achieved with instrumentation closer to the opening. 


ee i a a . 
betoibssq Liew ' est ai % 
Sia) om 

bna betoiberq —s t 

Yo siuess 6 yldsdor ‘ai 2ee361 nist 
ssade oidipilest .tebom ef to enoisit 

bao aiziit ni tluges qiebnued argeey 
eoinueasm oid anoyed yisdpite bes of 90 . 
q20d4 weled ot6 e9287 Gisiss e3o392 
eos? Lennes ed BA my them A 
eldipiipen eat nisps ita iriw prebaved’ pedte Co 

ed Biyaw oefdd yibiqes S70" sonemicad Deot eds 3 
bseoxe @9¢67 nietse bas on bile fem yptless n 
besewonow ftom sd Bluow vise) is-etedT . 

aaw 2if3 a6 doe stofvedea aot PsVA059 isi: 

.anoi ta7aqe eninge isi Hci teusake Atod pad" 
(sos beberle) arsitend ad? wodd. OF 2 bas ar.a vies 
priwollol bns one sub esis iievie {sitter § 


. 7 _ 


ed? .°-9M ahqwda ng deaqut tecepetb uo Ob ada 18 oxo 
iiew vidsnoesss ef. oon ‘@ineta tnibex 6&9 20 bn973 2 sa | 
eedar sid etd? of 10038 eatuod. 5.0 20 amid & aie a — 
si jostle eansvbs ssi sH7 s1oled 26 betamise: - 7 

ee 7 


«notiqut sini se 
fond gras tox muine Basnibena oto aoiv notsaterton 
, Prema eae be hermit 
smear a 7 ray @ | | 
eo 


’ 


139 


10? 


(°/./ hr) 


s} 
w 


\ 3 
1 ‘ MA 
We As 
Lane 

\ 

\ 


RADIAL STRAIN RATE 


io 
i 


10° 10' 10? 10° 
TIME (hrs) 


S 
n 


10° 


Figure 5.18 Comparison of Observed ana Predicted Radial 
Strain Rates - Test MC-7.08, Parallel to Principal Stress 
Directions 
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Figure 5.19 Comparison of Observed and Predicted Radial 
Strain Rates - Test MC-7.08, 45° to Principal Stress 
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Figures 5.20 to 5.23 present the predicted (shaded area 
or heavy line) and observed radial strain rates during and 
following the two tunnel widenings performed on Sample MC-7 
(Tests MC-7.10 and MC-7.13). Excellent agreement between 
prediction and observation was achieved. The divergence of 
the predicted and observed responses with increasing time 
can be related to the inaccuracy of the instruments at low 
Strains and the deviation, at that instrument location, from 
the mean time-dependent parameters employed. Once again the 
dominance of the face advance effect is evident. 

In both widenings extensometer #6 (6=135°) exhibits 
strain rates considerably in excess of the others and of 
that predicted. This excess rate resulted from the 
initiation of yielding or softening during excavation which 
occurred during the tests at this location (see Figures 
10.23 andi10.29) This yielding would have caused an upward 
shift (due to dilation and radial yield zone expansion) on 
the extensional strain rate curves resulting in the observed 
behavior. In the vicinity of the yielding zone the stress 
distribution is no longer elastic, as was assumed. The drop 
in radial stress will be greater than predicted and the 
Stress level is definitely greater than predicted (equal to 
1 if yielding occurred). Both these effects will results in 
an increase in the extensional strain rate. This increase 
may, in some cases, be sufficient to mask the face advance 
effect. A trend in that direction is evidenced by the 


response of extensometer #6 during the second widening 


‘ ey - 


ani 2 eniesarent Preapeanns ore 
vol ae osnaauxdani add to yossenosal edd 03 apie 


mort jaozduved snemetzent Ian2 34 somueacaliat sie mie: 


eddy cises son ,beyoigne aa anla J : 
y .tnebive @i joeiis ss - ay en So: 
atidtdas (*afi=}) a tesomoan®: 2: ia ae g 

Io bas etsiso etd Jo azsoxe enon! ae 
‘ad mo23 Begivest 4q3 easoKe eid? .bessibexq Lee -= 

doidw sotdsvsous pniaub gaineiioe to pribisiy to ne <— E 
scale eae ese) npttevol: its 35 s2a04 edd Pam da bs 

bi 6wou ne Beeuss ivi binow hethikatie aidT (es.ar 5 £S.t 

no (nolaasgze. snes Btity Letbe: bas noiislib ‘as | | 


bevisaedo sid mt entetongs gare etet nigra i 10. ' == 
way a 
gasts@ sad snot pn gd7 io ytinioty od: isd 4 a 


goth off .bemuezs sew ae rer regnel on @: sabia, 
ods Bas Bagnibeng madd zo450%— od {liw eemrze Teipeay 


| bee = 


a 
#} b 7 3 
: 
Ls rnb 7 
7 7 i 
7 7 . yy 
e = Li 
a me 


sonavon » exe) 4 enn 08 ausioitive sd 29869, 9 
akpwentb: asia, ai on han #3 


ah 294 _— at 


ene 


, a4 


140 


(°%./hr) 


S 


RADIAL STRAIN RATE 
S 


} 
=~ 


10:7 


10° 10° 10’ 10? 10° 


TIME (hrs) 


10" 


Figure 5.20 Comparison of Observed and Predicted Radial 
Strain Rates - Test MC-7.10, Parallel to Principal® Stress 
Directions 
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Figure 5.21 Comparison of Observed and Predicted Radial 
Strain Rates - Test MC>/. 107 845° to Prifieipalgstress 
Directions 
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Figure 5.22 Comparison of Observed and Predicted Radial 
Strain Rates - Test MC-7.13, Parallel to Principal Stress 
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Figure 5.23 Comparison of Observed and Predicted Radial 
Strain Rates - Test MC-7.13, 45° to Principal Stress 
Directions 
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(Fagure o, 23). 

According to Panet and Guenot (1982) the convergence 
rate depends on the extent of the plastic or yield zone and 
the ultimate tunnel convergence. As it has clearly been 
demonstrated that the face advance generally dominates the 
strain rate for a sufficient time after advance, it should 
be possible to use these measured rates to determine the 
extent of the yield zone. 

The tunnel convergence rate, as a function of the 
distance from the tunnel face, x, is given by Panet and 


Guenot (1982) as 


Unmex Pa=1G-2R* whe) Zk) ACR ZK Rte)? Eqntbs.12 
where R = radius of the plastic (yield) zone; 
k = 0.84 as before; and 


u(~)= ultimate tunnel wall displacement. 
According to Kaiser (1980) the ultimate tunnel wall 


displacement in an elastic, brittle plastic material is 


given by 
(haces 
Ute \r= PA RACT +a ie (R/a) -a="ipoa/ 2G Helis Sth 
e 
where r =a, (mein Lm to Zoo hl} 
e e@ 
a = dilation coefficient; 
m = tan acaomet ofa). 


@ = angle of internal friction; and 
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oO = unconfined compressive strength. 
c 
Rewriting Equation 5.12 in terms of time, t, and 


SUDSt1TCUL ING! dl t fore xen eLdss 
wilt) = [<2R2df wWSy7 ka k/CkRedea) | 2 Bans 5. 14 


where d' = rate of tunnel advance. 
Thee ratao! of? radial strain’ inthe’ elastae zonettor the 
displacement of the tunnel wall is given by: 


l+a 


e 
Gla) /ula)e= -F(R/r eta l/al2t(R/a) eke Tb sokebihs, Hela PS 
r 


and the ratio of radial strain in the plastic zone to the 
displacement of the tunnel wall by: 


+a 


etp 1 
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Ita 
teal 2tR/a) +a-1]} Eqns S816 


The above relationships allow determination of the radial 
Strain rate due to face advance in a yielding medium as: 


in the elastic zone: 


ed (tjt= (2R4(14+a)d' tle) 7krzal fk/(kRt+d' Covey 
G 


iartas 
[2(R/a) +a-1] EGn own 


ees way ws 
lige, Eee: ee 
of mp 

ene aesd ea ah 
eo nentienk 


wn in ae ae | 


Sf .ie 


ods od. aes sidney: wo so a ns 
nee Rimes igang? ef 39 


) Sua ts : 


— 


146 
in the plastic or yield zone: 


1+a 
Ctomantlaita)/2)+al(R/r) -1]}[2R2d'u(@) /ka] 
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Assuming that the plastic flow causes no change in volume 


(i.e. , @=)..0) sBquations™5.17 and 5.18 can be reduced to: 
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respectively. 

As mentioned previously, yielding was noted in the 135° 
orientation during widening of the tunnels of Sample MC-7. 
In an attempt to determine the extent of this yield zone, 
Equation 5.19 was employed to predict the relationship 
between strain rate and time for different extents of the 
yield zone. The actual behaviour cannot be modelled as 
Equation 5.19 assumes axisymmetric behaviour; hence, only an 
equivalent, annular, yield zone radius is reported. The 
results of this investigation are presented in Figures 5.24 
to 5.27 for the two extensometers located in this 


orientation. While these two instruments were located in the 
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Figure 5.24 Comparison of Radial Strain Rates: Observed and 
Predicted Face Advance - Test MC-7.10 (r/a=1.77,0=135°) 
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Figure 5.25 Comparison of Radial Strain Rates: Observed and 
Predicted Face Advance - Test MC-7.10 (r/a=2.12,0=135°) 
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Figure 5.26 Comparison of Radial Strain Rates: Observed and 
Predicted Face Advance - Test MC-7. 13 (r/a=1.36,0=135°) 
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Figure 5.27 Comparison of Radial Strain Rates: Observed and 
Predicted Face Advance - Test MC-7.13 (r/a=1.61,6=135°) 
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yield zone Equation 5.19 was still employed due to the 

uncertainty in the value of the dilation coefficient and the 

fact that the yield zone was localized (i.e., not annular) 
and therefore more rigorous analysis would not be justified. 

The results suggest that an equivalent radial yield 
zone would extend some 1 to 2 tunnel radii from the tunnel 
wall. The extensometer closest to the tunnel wall 

(Figures 5.24 and 5.26) predicts a smaller extent of the 

yield zone, however, this is probably due to the fact that 

there is less dilation occurring in this region as compared 
to the active yield front. In any event, the results of this 
study show that it should be possible to determine the 
probable extent of the yield zone around an advancing tunnel 
by means of displacement rate measurements. Limitations to 
be considered in the analysis are the certainty of the 
values selected for the variables in the equations and the 
axisymmetric nature of the solution. 

Several conclusions based on the above results may be 
drawn: 

1. Reliable prediction of material behaviour can be 
achieved by relatively simple models in the absense of 
yielding or softening; 

2. In the absence of yielding, immediately after the tunnel 
face advances past a given point, the radial strain rate 
is dominated by the continued advance, for a period of 
time that depends on the material properties and the 


rate of advance; 
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3. The face advance effect may not be as dominant in the 
presence of a propagating yield zone but, as will be 
seen later, will dominate for a longer period of time; 

4. Significant deviation from anticipated behavior is 
indicative of yielding or softening processes; and 

5. A first approximation of the extent of the yield zone 
can be achieved by analysis of the strain or closure 
rates due to face advance. 

In order to more fully understand the significance of 
the face advance effect the behavior under field conditions 
for this rock type has been investigated by means of the 
previously described relationships. The radial strain rate 
was determined at a distance of 1/2 tunnel radii from the 
tunnel wall (r/a=1.5): the tunnel was assumed to be 5 m in 
diameter and the field stress to be 12.5 MPa (corresponding 
to a depth of 500-600 m). It was further assumed that the 
material behaved in a linear elastic manner, the tunnel 
excavation process was not interrupted nor terminated before 
1000 hours and there was no scale effect on the material 
properties. f 

The radial strain rates experienced in the rock mass 
upon excavation of a tunnel at three rates,,.0.1, 1.0 and 
10 m/hr, were predicted and are presented in Figure 5.28. 
Both the total strain rate and the rate excluding the effect 
of face advance are plotted. The dominance of the face 
advance effect over a significant time interval (and hence, 


tunnel length) is evident. These results show that the 
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Figure 5.28 Comparison of Predicted Radial Strain Rates 
Assuming Various Tunnel Advance Rates in Wabamun Coal 
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effect of face advance on the measured strain rate cannot be 
neglected in this material until the face is further than 10 
tunnel diameters from the point of measurement. Displacement 
rate meaSurements made to evaluate tunnel performance or 
rock mass properties (for extrapolation) will therefore 
require adjustment prior to comparison with prescribed 
limits until the face is far advanced. In addition the use 
of the slope of the strain rate decay curve in estimating 
the long term displacement will result in significant 
underestimates of closure if taken prior to the dissipation 
of the face effect. 

Bigureso. 29epresetmussthemaccumulated strains versuse Nog 
time plots resulting from the excavations described above. 
From this figure it 1S apparent that the radial strain 
actually measured in the field will be very much dependent 
on the rate of tunnel advance as vel the time elapsed 
prior to instrument installation. If the instrument was 
installed 1 hour after the face advanced (Pt. A) 59 % of the 
total time-dependent rock mass response would be recorded 
regardless of the rate of advance. The amount of Strain 
recorded due to the face advance would be very much 
dependent on the rate of advance with only 3 % of the total 
aimaneadvance rate of M0em/noy45.98% at a. Uem/7hrgand, 91.1 % 
at 0.1 m/hr. If, however, the instrument was installed 
5 hours after the face advanced (Pt. B) the time-dependent 
rock mass response measured would be only 46.1 % of the 


total and the face effect strain would be 0.2 % of the total 
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atganeadvancesratesofusOem/hrasarS Yaa 400 mZhbrneand.65.26% 
at 0.1 m/hr. These results demonstrate the combined effect 
of instrument installation time and rate of tunnel advance 
on radial strain measurements and hence, instrumentation 
design. 

If the purpose of the instrumentation is to provide 
information regarding the possible extent of yielding and 
the development of radial strain around the opening then it 
should be installed very soon after face advance and be of 
sufficient range to accomodate large strains. If, however, 
the purpose is to determine the long-term strain to give an 
idea of ultimate lining load then more accuracy is required 
and late installation is acceptable. This generally implies 
a shortened range and therefore the instrument should 
possibly not be installed until the face advance effect has 
dissipated. | 

The degree to which the face advance effect influences 
the instrumentation design and analysis depends very much on 
the material through which the tunnel is excavated. As the 
face advance effect is directly proportional to the ultimate 
tunnel closure (excluding creep closure), the greater this 
closure the greater the influence of the face advance on 
measured strain or closure rates. Large values of the 
ultimate tunnel convergence may be a result of soft rock 
(low deformation modulus) or the creation of a yield zone. 
The influence of the yield zone on the face advance effect 


for the situation previously described is presented in 
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Figure 5.30 for the tunnel closure rate and Figure 5.31 for 
the radial strain rate. An increase in the radius of the 
plastic zone causes an increase in both the magnitude of the 
deformation rate and an extension of the time period over 
which it is of influence. This is clearly seen in 

Figure 5.32. In the elastic case (R=a) convergence to the 
ultimate is achieved within approximately 2 tunnel 
diameters. For increased extents of the plastic zone, the 
total closure increases as well as the distance from the 
face where it reaches its ultimate value. 

The significance of the face advance effect on the 
meaSured radial strain or tunnel closure rate has been 
clearly demonstrated and the expressions for its 
determination verified in the model tests. As these tests 
were of limited time and extent, extrapolations to field 
behaviour were made to provide additional insight. The 
results of this study have demonstrated that the influence 
of the face advance effect is very much dependent on the 
material through which the excavation is made and the rate 
of advance. Instrumentation for tunnel monitoring requires 


an appreciation of this effect for its reliable design. 


5.3.2 Development of Radial Strains Near the Tunnel Face 

The layout of the instrumentation provided an excellent 
opportunity to investigate the near face behavior during the 
initial excavation and subsequent widenings of the tunnel. 


The tunnel excavation and widening of Sample MC-6 were 
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Figure 5.30 Variation in Tunnel Closure Rate with Radius of 
Yield Zone 
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Figure 5.31 Variation of Radial Strain Rate with Radius of 
Yield Zone 
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performed in a relatively rapid manner whereas those 
undertaken on Sample MC-7 were performed at a slow rate, 
allowing more refined analysis of extensometer response. 
Sample MC-6 

Figures 5.33 and 5.34 show the development of radial 
Strain around the initial excavation of Sample MC-6. Little 
Straining was noted until the tunnel was within 1/2 diameter 
of the instrument level. Typically the instruments underwent 
a rapid increase in extensional strain as the tunnel 
advanced. Once the excavation extended beyond the instrument 
level, the slope of the radial strain versus time curves 
approached zero indicating little creep response (over a 
short time period). Instrument response was generally 
uniform with the exception of those located near the tunnel 
walt at) 6=90°,; [35° and 270° which exhibited approximately 
twice the average strain. This suggested that stress 
redistribution processes were occurring in these 
Orientations. 

Widening of the above tunnel to a diameter of 152 mm 
was done in a step-wise manner due to several interruptions 
in the drilling operation. The instrument response 
(Figures 5.35 and 5.36) generally followed the form of the 
tunnel advance curve with rapid increases in extensional 
strain during excavation. Upon interruption the straining 
abruptly levelled off into @ creep strain response. The 
advance of the tunnel face was tracked by the extensometer 


response. There was a wider variation in response as 
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Figure 5.33 Development of Radial Strain ~ 108mm Tunnel 
Excavation, Sample MC-6 (parallel to the principal stress 
directions) 
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Figure 5.34 Development of Radial Strain - 108mm Tunnel 
Excavation, Sample MC-6 (45° to the Principal stress 
directions) 
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Figure 5.36 Development of Radial Strain - Widening to 
152mm, Sample MC-6 (45° to the Principal stress directions) 
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compared to that noted during the original excavation, 
suggestinguthar stgnificant alteration of the sample shad 
occurred. Strains were maximum in similar, though more 
extensive, areas as during the initial excavation (45°-180° 
and-2/0%-375°) indicating that Stress? redistribution was 
continuing? 

Sample MC-7 

Figures 5.37 to 5.39 present the radial strain response 
recorded during excavation and widening of the tunnel for 
those extensometers nearest to and furthest from the tunnel 
wall. These instruments were always oriented in a true 
radial direction during all three tests. The complete record 
of measured radial strains for all extensometers during 
advance are presented in Appendix D. The time at which the 
tunnel face arrives at the level of each instrument is 
denoted by a vertical line, the type Be whe Spare ae to the 
appropriate radial strain versus time curve (also indicated 
by arrows). The near face behavior is also summarized in 
Tablets .3. 

In all tests, but particularly Tests MC-7.08 and 
MC-7.10, there was some increase in compressional radial 
strain prior to the expected extensional straining. This was 
caused by both an increase in radial stress ahead of the 
face (see Figure 5.14) and the time-dependent sample 
compression. 

Neglecting the compressional effect of face advance, 


the instruments near the tunnel wall (r/a<4.64) detected the 
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Table 5.3 Near Face Radial Strain Behavior 


Test No Location Face Advance | Face Advance | Max’m Radial Ea/€t 
Noted Complete Strain at (%) 


OPED Ce) SIS OS 


=) 5672 


A 
BSS 


-1.01+/-0.46 1.76+/-0.51 

e -0.47+/-0.13 2.32+/-0.65 
Sta.56to156 .16+/-0.82 
Sta.56to131 .38+/-0.62 


: : -O. : ; 
: : -0O. . : 
, : -O. y : 
3 : -O. : 2 
F 5 -O. 3 : 
: ; -O. : k 
: : -O. : . 
: : -O. : é 
: -0.54+/-0.16 1.25+/-0.52 
; SO) Siri SO). VE 0.90+/-0.30 
A 
Sta.56to156 a1 7+/—-O224 
Sta.56to0131 ; .13+/-0.23 


]=~NN=]NNN = 


-O. 
=O) 
=On 
=). 
BO 
=(6). 
Ole 
=(6)- 


a ee ee ee ee ee ae 
O9O0000000 


-0.67+/-0.23 0.60+/-0.29 51.7+/-18.2 
-0.58+/-0.12 ©, 45+/-0.02 61.5+/-2.3 
48 .5+/-13.8 
e 66.8+/-7.8 
Sta.56to156 .344+/-0.13 |55.4+/-14.7 
Sta.56to131 .344+/-0.14 |57.2+/-15.0 


= instruments near initial heading 
= instruments towards which tunnel expanded 
€a/et = ratio of strain ahead of the face to total strain 
NOTE: all values, unless otherwise specified, are given in 
terms of tunnel radii (a). 
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Figure 5.37 Development of Radial Strain - 42mm Tunnel 
Excavation, Sample MC-7 


169 


40 


face arrival ~ 


drilling interruption 


m | 
YL LLL 


( 


’ 


2.0 
TIME (hrs) 


-0.4 
-0.5 
-0.6 


N m 
° So 


(%) NIVYLS IVIOVY 


-0.1 


Figure 5.38 Development of Radial Strain - Widening to 
108mm, Sample MC-7 
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Figure 5.39 Development of Radial Strain - Widenin 
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presence of the initial tunnel advance approximately 1 
tunnel radius in advance while those further away (r/a26.4) 
detectedsitconly«0.5!tunnelsradii invadvancecon average: A 
Similar lag in response behind the advancing face was also 
noted; displacement levelled off at 1.76 radii for near 
instruments and at 2.32 tunnel radii for far instruments. 
This 1S indicative of the 3-D nature of the near face 
behavior. 

During widening, extensometers at lower r/a ratios 
(<2.62) reacted when the face was approximately 0.5 to 0.67 
radii from the measuring section. All instrument response 
levelled off between 0.45 and 1.25 tunnel radii behind the 
face. 

The distance over which the displacements occur due to 
face advance during the initial excavation is approximately 
Zatretunne lipradinyef or alt Petr ates aurea widenings this 
distance is between 1.59 and 1.18 radii. 

An increasing proportion of the total strain occurred 
ahead of the tunnel face during subsequent widenings. 
Approximately 20.6% (standard deviation s=17.5%) of the 
total strain occurred ahead of the face during the initial 
excavation, 133.5% (s=13.6%) eduning tthe pfurst lwidening and 
57.2% (s=15.0%) during the second and final widening (note: 
readings at Sta. 156 neglected because of proximity to 
sample boundary). The relative increase in the ratio of 
Strain ahead of the face during widening also was found to 


depend on the location of the extensometers. Those located 
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in the region towards which the excavation was expanded 
(denoted by 'e' in Table 5.3) showed increasing percentages 
of straining ahead of the face as compared to those near the 
imut taliveheadung (denotedebyd'n Phin) Table t5.3):. tThis 
resulted, most probably, from propagation of the yield zone 
in the direction of tunnel expansion. Appropriate values 
corresponding to the percentages of strain ahead of the face 


BGerauuvermto the stotal strain /for™the three tests are: 
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There was also a shift of the maximum strain rate from 
approximately 0.16 radii behind the face during initial 
excavation to 0.17 radii ahead of the face during first 
widening and then to 0.34 radii sheet! chee ne final widening. 
A probable explanation of the above monitored processes and 
their significance with regard to deformation monitoring 
Follow. 

During the initial excavation, the stress originally 
carried by the now excavated material is transferred around 
the excavation and to some extent to the rock at the 
advancing face. This increase in stress ahead of the face 
causes the initial increase in compressive radial stress and 
strain as previously mentioned. These increased stresses can 
lead to deformation and possible failure of the tunnel face 


resulting in the onset of extensional straining ahead of the 
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face. During tunnel widening, there is only a weakened (by 

initial excavation) unexcavated rock core (hollow cylinder) 

left to resist the additional stresses. Overstressing near 
the tunnel wall is increased. This results in the increased 
straining ahead of the tunnel face. This process as well as 
yielding of the tunnel walls ahead of the face has 
effectively advanced the tunnel ahead of its present 
location. 

The significance of the above findings to the 
instrumentaion designer are as follows: 

1. In order to maximize the amount of information on the 
deformation field around the excavation, where multiple 
heading advance is employed, the instrumentation should 
be installed from the pilot headings. With each 
successive expansion of the excavation the magnitude of 
the post-excavation deformation will decrease; 

2. Particular attention should be paid to deformations in 
the direction of cavity expansion. Yielding initiated 
around the original heading will tend to propogate in 
the direction of expansion as the magnitude of the 
stress change in the other direction will be much lower; 
and 

3. Interruptions in the tunnel advance provide an 
opportunity to determine tunnel stability by 
displacement rate criterion or to determine the 


time-dependent rock mass response. 
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6. CONCLUSIONS AND RECOMMENDATIONS 


6.1 Introduction 

A rational design of underground excavations in rock 
should employ the concept of the convergence-confinement 
method where the the ground-support interaction is 
considered as presented in Figure 6.1. The convergence curve 
(ground reaction curve) depends on the rock mass properties. 
Tt is linearsin a linear elastic material, curved if 
yielding occurs and may tend upwards if gravitational forces 
dominate in weak loosening rock. Equilibrium can be reached 
at any point along this curve if the appropriate pressure is 
applied at the corresponding displacement (intersection of 
confinement (support) and convergence curves). Where the 
rock mass properties are time-dependent the convergence 
CUunVemisha function Ofmtame, oveequatey Bivientt Mi eisy Kote baep an 
(u = 0, p = po) with time as shown in Figure 6.1. 
Equilibrium is therefore not assured at the intersection of 
the convergence and confinement curves at a given time 
(pt. A). Additional creep deformation can shift the 
equilibrium position to Pt. B on the confinement curve for 
infinite time. This may result in overstressing or failure 
of the support system. Similarly a shift of the convergence 
curve may occur if the field stress increases due to nearby 
mining operations. This may lead to a further increase in 
support load (Pt. C). For design it is therefore necessary 


to determine the ground convergence curve from the rock mass 


174 


ee Poa 


Ol ete bes 13 
‘ier rho llega’ 26, | 
Lvs os rotenone vena oan Seo 


sviup soaspyaveac> af? 148 stoi ni Awanatend ag bevebleties 

.esit7aqeaq e2em 1907 sid no zbaaqed fevivo noitossz-bnuemg) 

ti bavago detaaten sitaale ysenit 4 nk assnif ai 10 

asnio? ienoftepiveie 2¢ abrewqu’ S09 -yem Bae core 

bedoset 24 15D muitdiiiup® .i207 painezool jos at-adeninab 
ei swwete1wy eseiiqesqqe edt ti evaivo eid? pnole soho gs: wu : 
te nol¥oeeisdai) seemepslqeid paibaoqas220%- end ts bsifqas : 
od? ovadW «(eevee eoteptevno> brs (s10qGque) tnamentings 
eonepsevnos 3f7 gnebreqesb-enis ste: eettreqoxg 22am A203 ; 

aipixe aad tueda gnidstoz , emit 20 noisonu? Ss ei eviuo 

_tv@ evup#t ai avoda es emis driw (og = quo = a) 
to noitveeseint ad& 48 Betueae Jon a1otezeds ef auiuditivps 
onid nevig s Je eevmrs tnementines bas sonsp2evn0o odd 
od stile ng medaemr0i26 qevis LenoistbbA «(A 240), : 

soi evo snsmenttaco-edd co 4 .29 OF notsieoq avinditiupe 
eiulis? 10 grteestievevo of divest yom eidT .omis seiniiad 
sonepsevacn add 30 dtide/e ylaelimiz -mosaye sroque ats 20: 
ydtaen of sub esaseromi Resise Sisti sda 1 swore tam ovaUS. 
a szassoni. redaaui ® of peel yan sta enoisexego eninia 


RADIAL STRESS 


__original ground stress 


support curve 

linear elastic 
elasto-plastic | rock 
time-dependent 


stress change due to 
nearby mining operations 


A 


-~ loosening 


RADIAL DISPLACEMENT 


Figure 6.1 Confinement-Convergence Curve 
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properties and the established mode of opening behavior in 
order to calculate the shape of the confinement curve for 
the selected support system and support activation time. 

The ground convergence curve can be determined 
analytically by methods such as those presented by Brown et 
al. (1983) if the rock behavior is time-independent. For 
rock masses with time-dependent responses, the determination 
of the ground convergence curve becomes much more difficult. 
In this situation two approaches to the problem are 
possible; an analytical or an observational approach. 

The analytical approach entails employment of 
visco-elastic or visco-plastic material models. These, 
however, have difficulty in incorporating geologic 
complexities and different modes of failure. In addition, 
the required parameters can seldom be determined in a 
rational manner because they do not penemoete actual rock 
properties. Hence, this approach is generally unsuccessful 
in accurately representing the ground response. The method 
proposed by Ladanyi (1974) can, however, be used for 
Parametric studies. 

In the observational approach, the designer resigns 
himself to the fact that the ground response is analytically 
incalculable and that the design must be derived from 
in-situ experience. This experience is gained from the 
analysis of performance records of openings constructed 
under similar conditions or preferably during pilot 


excavations. 
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Deformation or displacements serve as the best 
indicator of opening performance because stress and stress 
changes are difficult to measure reliably. The displacements 
and especially displacement rates constitute the most easily 
obtained and most reliable measurement and, hence, are 
selected as performance indicators during most tunnelling 
operations. 

Because of the importance of the deformation monitoring 
program to the design of underground excavations in rock, a 
review of its application and limitations was in order. 
This, essentially, constituted the thrust of the research 


reported herein. 


6.2 Conclusions from Model Tests 

The design verification and stability assessment 
aspects of deformation monitoring in typical underground 
excavations are performed in different manners. The 
stability of the excavation is generally made on the basis 
of convergence measurements recorded at regular intervals 
along the tunnel axis while design verification is performed 
at selected sites representative of typical conditions by 
more intensive instrumentation arrays, particularly borehole 
extensometers. Both generally employ displacement rate 
criteria in their application, however, design verification 
also relies on displacement magnitudes (for serviceability 
requirements) and often includes support stress 


measurements. 
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The response of the selected instrumentation depends on 
the material property distribution and its stress level 
dependence. In a heterogeneous, discontinuous rock mass the 
variation in instrument response will be significant. The 
application of a universal displacement criterion in such 
Situations is not practical. Judgements must be based on the 
change in response at a given location or on the difference 
when compared to other instruments under similar conditions. 
Similar conclusions were also drawn by Guenot et a]. (1984). 

Convergence measurements were found to be insensitive 
to the deformation field around the excavation. The 
non-axisymmetric nature of the yield process was not always 
readily identifiable from the convergence records. In 
addition the magnitudes of the observed closures did not 
allow for the identification of the deformation process or 
its location and extent. pReeneoset eck however, were able 
to provide the required information concerning the 
displacement field. 

The magnitude of the deformation observed at any 
Station along the tunnel axis will depend on the time of 
instrument installation and the stage of the excavation 
process. In order to maximize the deformation information, 
instrumentation should be installed ahead of the advancing 
face. This is often impossible or impractical in deep 
excavations. In these instances, installation at the face 
will provide sufficient information for analysis. In 


multi-stage excavation, instrumentation should be placed 
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during the heading excavation as the magnitude of the 
deformations decrease with each successive stage (unless 
instability arises). 

The intended use of the instrumentation, to a large 
extent, dictates the measurement range and time of 
installation. If the objective is the determination of the 
Gevelopment of radial strain around the opening in order to 
ascertain the location and extent of yielding or the 
verification of design assumptions then the instruments 
should preferably be installed ahead of the face or as soon 
as possible after face advance and be of sufficient range to 
accomodate large strains. If, however, the purpose is the 
prediction of ultimate lining loads from extrapolated creep 
Strain curves then later installation of reduced range (more 
accurate) instrumentation is permissible. In fact, it may be 
appropriate in these circumstances Loqanhen dissipation of 
face advance effects prior to installation. 

The distance (in tunnel diameters) over which the 
displacements are influenced by the advancing face depends 
on the type of excavation as well. The distance is maximum 
during the initial excavation and becomes progressively less 
in subsequent widenings. In such circumstances, instruments 
designed to detect the yield zone must be installed during 
the initial heading excavation or at the face during 
widening. Instrumentation for long term displacement 


measurements may be installed soon after final widening. 
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In order that the maximum amount of information be 
obtained from the instrumentation frequent readings are 
required. This is becoming more practical with the advent of 
rugged, remote data acquisition systems and computer data 
processing. Continuous records of displacement can yield 
information on the intermittent, localized behaviour which 
may be of assistance in understanding the overall tunnel 
behaviour. Frequent readings near the tunnel face are 
essential if determination of the face advance effect is 
required. 

During external loading tests, it was observed that 
critical deformation rates, indicative of a behaviour mode, 
existed for the small tunnels in the coal. These critical 
rates, which were employed as tunnel behavior evaluation 
criteria, were applicable to the time-dependent rock mass 
deformation rate (i.e., creep rate). The existence of such 
Critical rates should be tested by field measurements. 

The deformation rate measured during tunneling 
operations is a function of the rate of tunnel advance and 
the time-dependent rock mass properties. Time-dependent 
measurements made to evaluate tunnel performance or rock 
mass properties therefore require adjustment prior to 
comparison with prescribed limits until the face is far 
advanced. In order to apply the critical rate functions 
previously described, it 1s necessary to separate the 
observed rate into that due to the face advance and that due 


to creep including time-dependent yield processes. This was 
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achieved by extension of an expression relating tunnel 
closure to face advance proposed by Panet and Guenot (1982). 
Relationships for both the tunnel closure and radial strain 
rates due to the advance of the tunnel as functions of time 
and advance rate were derived and subsequently verified in 
the model tests. Substitution of reasonable values for the 
parameters in these expressions demonstrated that the face 
advance may dominate the observed deformation rate for a 
considerable distance behind the tunnel face (approximately 
10 tunnel diameters for this coal). Where the slope of the 
Strain rate decay curve is used in estimating the long term 
displacement, elimination of the face advance effect must be 
achieved in order to prevent signifigant underestimates of 
tunnel closure. 

Closure rates in excess of the closure rate functions, 
while indicating a behavior mode, did not indicate where the 
process was occurring. Measurements of radial strains were 
found to be necessary in combination with closure rates to 
delineate the location and extent of the yielding or 
softening zone. Furthermore it was found necessary to 
measure the tunnel closure rate or radial strain rate in 
many directions and continuously, beginning as soon after 
the tunnel advance as possible, to ensure correct 


interpretation of the mode of abnormal behavior. 
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6.3 Recommendations 

Deformation monitoring of underground excavations in 
weak rock can be adequately achieved by the use of 
convergence meaSurements and borehole extensometers. These 
two means complement each other and should be incorporated 
as such in the monitoring scheme. This research has shown 
that the onset on incipient instability can be predicted by 
comparison of tunnel convergence rates with appropriate 
critical rate functions but that the location and extent of 
the instability requires the use of borehole extensometers. 

A recommended procedure for the suitable application of such 

a deformation monitoring scheme follows: 

1.. Establish critical closure rate functions. These may be 
obtained from performance records of tunnels of similar 
geometry, constructed in similar geologic settings or 
from test sections or pilot headings of the same tunnel; 

2. Take convergence readings as soon as possible after 
tunnel advance. Frequent readings should be taken until 
the mode of behavior is well established. 
Multi-directional closure measurements should be made 
with a 4-point convergence array as a minimum; 

3. Separate the time-dependent closure from the face 
advance effect by application of the derived expressions 
and compare with the appropriate functions to determine 
behavior mode; and 

4, I£ yielding is indicated or anticipated, install 


multi-point borehole extensometers around the excavation 


ae Se 
rm | _ “a 1 eal ve es, 
» j hate irs ite 


bean Qake yeti bleu ; i 7 — saad 
at enoiseveoxs 5 Kstinem cobs te 


peep 5) oe 


“is rane ms 
saedt crcasascenssi® atofored f eT, —_ 
petaroqrosnt ed bluode bane 29d30 inenplanee: i one 


nveda. ead, dgresger ai dt saanionset inc 6 —_ 

ed Rassiderg ed,nge tbbidatent sasigtont ao edacaade 

agotaqaxaqe dam aods2 scaspyevndo tenqu? 20 npatasqmos 

14 ghedxd be soisezel 2d3 teddy Jud, anoksany 938% Kapkaine’ 

.eisismagngsxe Ssiodaied Yo pau os) apa eesiupoy ysbtidesants edd 
qous 2% nottesilqqa eldasiue stg, 10%, ssubeccng babnommoaenek 
_rawol{o2 america. pmizosinon nottemtoleb a 

ed, yam sasd? .enoitonus eter e2u2e0lo fepiiias deiidssem tad 

ssiimte to alennuy to abtoses somsmsotieg oor. benkeddon a) 

ro apnidiee sipelobp zakimte ni) beavinsenoo, weizemosg i 

:Lanaud smpe aria Bo. .apadbeort toliq 2© enoti 2e2: a 

ie43n, sidieeog ea acoe as epnibse2 eoneprevnom 

[itau nodes od Blues) eppribees, Jneupes4 ,eonsvba:f 

. , sbedetidetes Lfow 2i z0lveded to 

sbhsm od Sluoda sipegasuasen, swaolp Lenotsoesib= ha hums 


- 


183 


or locally, if the probable location of yielding is 
known, to delineate the location and extent of the 
overstressed rock so that remedial measures may be 
taken. 
The above procedures are very much dependent on quick data 
assembly and processing as well as flexible contract 
arrangements that allow for instrumentation based on 
rational design criteria rather than predetermined criteria. 
However, instruments may be placed for other reasons as 
mentioned in the introductory chapters. 

This proposed deformation monitoring program was 
derived from a series of model tests conducted on coal 
Samples under controlled laboratory conditions. To verify 
their suitability for other materials would require 
additional laboratory investigations. However, it would be 
much more beneficial to explore their validity in field 
scale operations and possibly some parametric studies 
performed with sophisticated numerical material models. 

An additional requirement identified in this research 
is the need for the development of instrumentation capable 
of continuously recording the radial displacement or strain 
ahead of the face in deep excavations. Only with these fully 
defined deformation records will verification of design 


assumptions be established. 
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8. APPENDIX A - DETERMINATION OF TIME-DEPENDENT PROPERTIES 

Linear Elastic, Isotropic Parameters 

The Young's modulus, E, and Poisson's ratio, JU, 
required for the linear elastic, isotropic representation of 
material response were determined by one of the following 
methods depending upon the sample loading history. 
Method 7 

In this method, developed by Guenot (1979), the 
response of an extensometer during the application of an 
increment of stress before and after tunnel excavation is 
compared. It assumes that, during subsequent reloadings of 
the sample, the elastic parameters are loading history 
independent, the rate of loading is similar for both 
situations and that linear elastic solutions and 
superposition are applicable. 

The radial strain increment (Ae ) due to a change in 
ehesma jor principal) faeld stress ie inechevintact 


Vv 
specimen under plane strain conditions is given by: 


Ne ‘= ‘AcinGieV Ph Ghe2U) Cian) AOI-Niicos26 )/2E 8.1 
r Vv 
where N = ratio of minor to major principal field stress 
and 
6 = angle with respect too . 
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Equation 8.1 can be rewritten to yield an expression for 


Young's modulus as: 
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E = Ao (1+¥V)[(1-2U)(14+N)+(1-N)cos20]/2Ae ae 
Vv is 


Both Ao and Ae are measured during testing and can be 

Vv u 
incorporated into a single variable 'm' which represents the 
slope of the tangent to the stress strain curve at a given 


Sie CSS. 


m= Ao /e eae 
v is 


(eS) 


Equation 8.2 can then be rewritten as: 
E = 0.5 m (1+V)[(1-2V) (1+N)+(1-N)cos26] 8.4 


Under conditions of equal lateral field stress (N=1), 


Equation 8.4 can be rewritten as: 
BE Som(aev) (4=2V) SR 


The radial strain increment due to a change in the major 
principal field stress in a specimen containing a tunnel of 
radius 'a' under plane strain and equal lateral stress (N=1) 


conditions is given by: 


> 
ia) 
tN 


Ao [(1-V) (1-a2/r?2)-U(1+V) (1+a?/r?)]/E 8.6 
r Vv 


where r distance from tunnel centre to point of 
measurement. 


Equation 8.6 can be simplified and rewritten to yield an 
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expression for Young's modulus as: 
EB = m(1+V)(1-2U-a?/r?) Ba 


Substituting m, for mein Bquation 6 .57anG mM. in Equation. s.7 
and then equating the two expressions yields an expression 


for Poisson's ratio as: 
V = [m, (1—-a? /r?-my ) ]/ 24mm, ) 8.8 


The value thus obtained can be substituted in either 
Equation 8.5 or 8.7 to yield the value of Young’s modulus. 
Method 2 

In this method comparison of the response of an 
extensometer under different loading conditions (N values) 
on the intact specimen are compared. Conducting two tests at 
different stress ratios, N; and N2, and employing 
Equation 8.4 results in the following expression for 


POUSSOn S racio: 


V= 0.5 - [m2(1-N2)-m,(1-N,) ]cos2é6 8.9 


rT 


2{m,(1+N,)-m,(14+N2)] 


This value may then be substituted into Equation 8.4 to 


determine Young's modulus. 
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Linear Elastic, Transverse Isotropic Parameters 
For a linear elastic, transverse isotropic body, whose 
plane of isotropy is the r-@ plane the relationship between 


stress and strain is: 


( WE, Yrg/Ep-Uzp/Ez 0 0 0 Ge 
en psa Wis OL ase VS 0 0 Ge 
Cy \ | er/Er Yer/Er VEz 0 0 0 Gs an 
fos 0 0 0 Ver © 0 fhe a 
yee 0 0 0 OMNI /2G>,00 0 Te 
yee 0 0 0 0 O  WUre/Er | | 7 zr 


where e denotes the normal strains, y the shear strains, o 
the normal stresses and 7 the shear stresses. The 5 


constants required for material description are: 


E = the Young's modulus in the plane of isotropy Ges) 
5g 
plane); 
E = the Young's modulus normal to the plane of 
Zz 


isotropy (z direction); 

G = the shear modulus relating shear stresses to shear 
Strains out of the plane of isotropy; and 

v,vV = the Poisson's ratio relating straining in the 
direction of the second subscript to the strain in 


the direction of the first. 


Poisson's ratio and Young's modulus in and normal to the 


plane of isotropy are related by the following expression: 
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Considering normal stresses and strains only and 


Substituting Ko for o the following expressions are 


@ 6 
derived: 

eel=n[oe(tekVi in] VerckB)/E Gee 2 2. 
ig r ré feed. Pe r 

e = [o (K-U ) - UV o J/E Glos 
6 G ré aya Ae r 

e '= [a —-V of €1+kK)]/E 8.14 
Z 2 ae rae Z 


where: K = [(1+N)-(1-N)cos26]/[(1+N)+(1-N)cos26]. 


Under plain strain conditions (e = 0) Equation 8.14 can be 
Zi i 
rewritten to yield an expression for U as: 
zr 
UPe ioe lor Gieky) he ee 
zr Zz r 
Under conditions of hydrostatic compression (0 = 90 =@ = 
8 6 


ao ) Equation 8.15 can be rewritten to yield an expression 
Zz 
ECOren Once .U is known as: 
Zz vale 


BE = 0(1-2U )/e eke Se 
Zz yA 9 ied 


Equation 8.12 can then be rewritten as: 
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tests at different stress ratios (N) will result in two 


equivalent expressions for E . This then allows for the 
r 
determination of VU by inputting the values of the measured 


20 
variables in the following expression: 
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Once VU is known, substitution of its value in 


ré 
Equation 8.17 yields the value of E . Substitution of the 
ie 
Valuesmior Ul shmeandel intosequationes el iay2eldsethe 
ZG Zz r 
value of UV. 
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Similar expressions for the relationships and constants 


can also be derived in cartesian coordinates. 
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9. APPENDIX B - RADIAL STRAIN AND TUNNEL CLOSURE RATES 
MEASURED DURING CREEP TESTS PERFORMED ON THREE SAMPLES 
The mean and range of the radial strain rates measured 

by extensometers at times of 1 (empty symbol), 5 (half full) 
and=20-hours (full symbols) aftter=toad~applicationare 
presented in Figure 9.1 for Sample MC-5, Figure 9.7 for 
Sample MC-6 and Figures 9.12 to 9.14 for Sample MC.7. Also 
shown on these figures are the boundary strain rates defined 
as the ratio of sample loading boundary displacement to 
sample width per unit time. The remaining figures present 
the measured tunnel convergence (closure) rates; Figures 9.2 
tomoromrlor Samp] cm MG—»> . rrgures 9.0 tO oot for MC-6vand 
Figures 9.15 to 9.19 for Sample MC-7. Also shown on these 
figures are the two behavior evaluation funct1o0ns, CCR and 
PCR (range denoted by shaded area), at a time of 20 hours 
afterload application.=The region denoted by the arrow 
between the CCR and PCR denotes the range of closure rates 
where behaviour is determined by the trend of closure rate 
versus time curves (in a double log plot). If the closure 
rate function converges to the CCR, yielding has terminated 
and a new stable equilibrium will be reached. If the rates 
are approaching the PCR range then yielding will tend to 
propagate. A stable equilibrium may still be reached but the 


accumulated closure may be unacceptable. 
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Figure 9.1 Strain Rates: MC-5.02 and 5.03 
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Figure 9.2 Closure Rates: MC-5.05 
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Figure 9.3 Closure Rates: MC-5.07 
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Figure 9.4 Closure Rates: MC-5.09 
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Figure 9.6 Closure Rates: MC-5.13 
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Figure 9.7 Strain Rates: MC6.01S and 6.02 
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Figure 9.8 Closure Rates: MC-6.02, 6.03, 6.04 and 6.05 
(108mm tunnel) 


* 


wl wahacarnny 


’ 
a 
- 


, 4 
- — en ee dae oe ey a) F 
— SS, —, 6k ps = + tae . 
— | ; : Palle 
r nn ances - maui see err 4 
: . ? a) mien a. —— ) 
3 ae — 2 th Z ; ath he 
, i = — oe — -_ ook ~6= 6 orm & ——— > 
f : > Wiggs som, 2 : i 
ony j 
re _. es = a ee er ee od 
on a { ee ee eel > 
s I. * -_ s 4 7] 
- ————- om , ; 
ae 3 ane ree 
x Pod 
: »— ——- ~~ ¥ 
— a ee es — nay bse aaa — fo 
: a «aa . ie * é oS -@ (he es Pus 
es --- < > @ oases; a 
t — oot 13 Ee 
— a ee aes ? 


ap a tmanre 
‘ 
! 
t 


ee ub te 


<3 Oe oe 


a eae Kagawa 4 ie . - aoe ae “ ; 

= Gis =o Le Senne +4 
-— >? are oe eT ae “—— 4a" 
2 “ # 1 +- ada - Sid > ~ = 7 


beeen aee Comes! ie 
i - : es 


eS 2) oo — | 
a ee na i FC Pe i yaa nea: io aR tie 
—s = < = }~ ttn pl r ee — ~-* whe oy o> 7 


— 
Daten mec 
a <8 | Oe = a a as 


’ 
eS t ——— to ——— ——- a 


“a. 47 * _—T  s 
S.* - pam -, al 


eS i ale 
stm a= SS. stare) eee we pdm o< 
a 


a 


7 


205 


if 


Figure 9.9 Closure Rates: MC-6.05 (152mm tunnel), 6.06 and 
6.08 
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Figure 9.10 Closure Rates: MC-6.09 
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Figure 9.11 Closure Rates: MG=6..10 fand: 6.4718! 
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Figure 9.12 Strain Rates: MC-7.02, 7. 
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Figure 9.13 Strain Rates: MC-7.06 
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Figure 9.14 Strain Rates: MC-7.07 
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Figure 9.15 Closure Rates: MC-7.10, 7.11 and 7.13 
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Figure 9.16 Closure Rates: MC-7.12 
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Figure 9.17 Closure Rates: MC-7.15 
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Figure 9.18 Closure Rates: MG—vi in, 


214 


4 


7. #% stetrs, 


ei = — ag Se 
ae, he 


1 


" . 


‘bel ate ye 


meh ae 


aivact @ 
maa 6 


Ly 


* 


ied 
—— 


~ - 


215 


4d Neoge | - bE 


4 
4 


| MC 


Figure 9.19 Closure Rates: MC=7.037 07. (¢9anda7 216 


: eet bee 
ee ae 


10. APPENDIX C - TUNNEL CLOSURE AND RADIAL CREEP STRAIN 
INCREMENT PATTERNS OBSERVED DURING TESTING OF THE THREE 
SAMPLES 

The time-dependent tunnel closure (across 4 diameters) 
and radial strain observed during a single creep stage 
(extensional strains and outward movement of the tunnel wall 
shown negative) are presented in Figures 10.1 to 10.6 for 
Sample MC-5, Figures 10.7 to 10.22 for Sample MC-6 and 
Figures 10.23 to 10.39 for Sample MC-7. The contours of 
radial strain were generated by a computer program that 
neglects the tunnel boundary and hence, do not accurately 
reflect strains near the tunnel wall. They do however 
adequately represent strain contours in the Vicinaty Of she 
extensometers. Only those stages where significant creep 
strains were recorded are presented. These figures assist in 
the determination of zones of active overstraining and 


possible overstressing. 
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Figure 10.1 Creep Closure and Contours of Creep Strain h- 
MC-5.05 at 12.8 MPa 
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Figure 10.2 Creep Closure and Contours of Creep Strain - 
Mc-S.07) ate 12,29 Mra 
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Figure 10.3 Creep Closure and Contours of Creep Strain - 
MC-5.09 at 13.4 MPa 
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Figure 10.4 Creep Closure and Contours of Creep Strain - 
MG=53 14 at let MPa 


727 


Figure 10.5 Creep Closure and Contours of Creep Strain - 
MG-Salnhat 13.0 MPa 
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Figure 10.6 Creep Closure and Contours of Creep Strain - 
MC 25°. 13) at wih. 0 MPa 
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Figure 10.7 Creep Closure and Contours of Creep Strain - 
MC-6.02 at 12.6 MPa (t=28.8) 
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Figure 10.8 Creep Closure and Contours of Creep Strain - 
MC=6- 02 ait, Oli2 46; MPa 4(b=29.9) 
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Figure 10.9 Creep Closure and Contours of Creep Strain - 
MC-6.03 at 12.4 MPa 
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Figure 10.10 Creep Closure and Contours of Creep Strain - 
MC-6.04 at 14.9 MPa 
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Figure 10.11 Creep Closure and Contours of Creep Strain - 
MC-6.05 (108mm) at 12.4 MPa 
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Figure 10.12 Creep Closure and Contours of Creep Strain - 
MC-6.05 (152mm) at 12.4 MPa 
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Figure 10.13 Creep Closure and Contours of Creep Strain - 
MC-6.06 at 12.4 MPa 
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Figure 10.14 Creep Closure and Contours of Creep Strain - 
MC-6.07 at 12.4 MPa 
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Figure 10.15 Creep Closure and Contours of Creep Strain - 
MC-6.07 at 14.8 MPa 
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Figure 10.16 Creep Closure and Contours of Creep Strain - 
MCG-6..08 at) T2a5' MPa 
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Figure 10.17 Creep Closure and Contours of Creep Strain - 
MG-6.:08 at, 15.0 MPa 
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Figure 10.18 Creep Closure and Contours of Creep Strain - 
MC-6.09 at 12.4 MPa 
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Figure 10.19 Creep Closure and Contours of Creep Strain - 
MG+6 09 gat 115 70 Mpa 
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Figure 10.20 Creep Closure and Contours of Creep Strain - 
MC-6.10 at 12.1 MPa 
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Figure 10.21 Creep Closure and Contours of Creep Strain - 
MC-6.14"abl 42.0 "MPa 
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Figure 10.22 Creep Closure and Contours of Creep Strain - 
ME-Guiieat fomosMrs 
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Figure 10.23 Creep Closure and Contours of Creep Strain - 
MC-7.10 at 12.5 MPa (t=31.5) 
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Figure 10.24 Creep Closure and Contours of Creep Strain - 
MC-7.d0cat 412.5 MPag(t=91.5) 
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Figure 10.25 Creep Closure and Contours of Creep Strain - 
MG=7 si tiat l2eo MPa 
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Figure 10.26 Creep Closure and Contours of Creep Strain - 
MG=7 Si2 gat }h2 35eMBRa 


243 


SIDE 4 


SIDE 1 
SIDE 3 


SIDE 2 


Figure 10.27 Creep Closure and Contours of Creep Strain - 
MC-7.12 at 14.3 MPa 
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Figure 10.28 Creep Closure and Contours of Creep Strain - 
MC-7.13 at 12.5 MPa (t=20.0) 
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Figure 10.29 Creep Closure and Contours of Creep Strain - 
MC=7.13 at 12.5 MPa (t229.7) 
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Figure 10.30 Creep Closure and Contours of Creep Strain - 
MC-7.13 at 12.8 MPa (t=89.7) 
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Figure 10.31 Creep Closure and Contours of Creep Strain - 
MC=7.1% at 1295eMPa 
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Figure 10.32 Creep Closure and Contours of Creep Strain - 
MC-7.15 at 72 59MRa 
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Figure 10.33 Creep Closure and Contours of Creep Strain - 
MC-7.15 at 10.0 MPa 
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Figure 10.34 Creep Closure and Contours of Creep Strain - 
MC-7.15 at 12.5 MPa 
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Figure 10.35 Creep Closure and Contours of Creep Strain - 
MC-7.15 at 14.0 MPa 
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Figure 10.36 Creep Closure and Contours of Creep Strain - 
MC=7.916 fat ‘lize SeMPa 
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Figure 10.37 Creep Closure and Contours of Creep Strain - 
MC—/ AVP at, 7AOtMPa 
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Figure 10.38 Creep Closure and Contours of Creep Strain - 
MC—7/.,li7) ate 10).,0) MPa 
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Figure 10.39 Creep Closure and Contours of Creep Strain - 
MC—/ elie Leelee OeMoa 


11. APPENDIX D - DEVELOPMENT OF RADIAL STRAINS NEAR THE 
TUNNEL FACE - SAMPLE MC-7 
The three sets of figures presented herein, show the 
radial strain development around the advancing tunnel face 
for the conditions of initial excavation (Test MC-7.08) and 
subsequent widenings (Tests MC-7.10 and 7.13). The initial 
figure presented in each set shows the position of the 
advancing face as a function of time. The time at which the 
face arrives at a particular instrumentation station along 
the tunnel axis is given by the intersection of the advance 
curve and the dashed line corresponding to the desired 
Station. The remaining figures in each set show the radial 
strains measured by the extensometers during each test. A 
plan view (to scale) of the extensometer positions is 
presented in the upper right corner of each figure. The time 
of face arrival and tunnel completion are noted on each 
figure (dashed lines). 
The excavation of the 42 mm diameter tunnel 

(Test MC-7.08, Figure 11.1) was interrupted due to 
mechanical problems at Stationen46s for atperiodgor 
0.3 hours. After this interruption the rate of advance was 
increased from the previous 0.026 m/hr to 0.102 m/hr (on 
average). Little extensional straining was observed until 
the tunnel was within 0.5 diameter of the instrumentation 
station. Typically the extensometers exhibited a rapid 
increase in extensional strain as the face advanced with a 


subsequent rapid decline. Virtually all straining was 
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complete within a diameter behind the face. This may be a 
result of the model's boundary conditions. Extensometers at 
Station 81 (6=135°), Station 106 (@=45°) and Station 156 
(@=45°) and Station 156 (6=90°) exhibited straining much in 
excess of others at the same distance from the tunnel wall, 
however, the magnitudes of the strains were still very small 
Oe). 

Widening of the excavation to a diameter of 108 mm 
(Test MC-7.10, Figure 11.9) was interrupted twice by 
mechancial problems. The first occurred at Station 146 for a 
period of 0.43 hours and the second at Station 170 fora 
period of 0.23 hours. The rate of advance, initially at 
0.071 m/hr was increased to 0.984 m/hr after the second 
interruption. Strain patterns were similar to those observed 
in Test MC-7.08, however, extensional straining was 
generally not observed until the face was within 0.26 
Giameters of the instrumentation station. The distance over 
which the greater part of straining occurred (in terms of 
tunnel diameters) was much reduced. Strains greatly in 
excess of the average were observed in the 45° orientation 
@ndeparticulanly in thes )35e orientation. 

Widening to a diameter of 152 mm (Test MC-7.13, Figure 
11.17) was performed without interruption at a rate of 0.073 
m/hr. Strain patterns were similar to those observed 
previously with continued higher than average straining in 
the 45° and 135° orientations. Some high strains were also 


observed in the 90° and 225° orientation close to the tunnel 
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Figure 11.1 Tunnel Advance Test MC 
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Figure 11.6 Development of Radial Strain - Test MC-7.08 
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Figure 11.9 Tunnel Advance - Test MC 
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Figure 11.17 Tunnel Advance - Test MC 
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Figure 11.18 Development of Radial Strain - Test MC-7.13 
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Figure 11.21 Development of Radial St 


mo 

™~ 

| 

O 

> 

a) 

149) 

oO 

EB 

! 

CG 

“a 

ice} 

» Se 
Cama 
o 
— 
—) 
ie) 
=) 
© 
Q 
~ 
\O 
<= 


Sta. 


i 
ye TT 


| 
i 
i 


0.€ 


oparecer see" 


2." 


i ; 


ee ma ll lt ee a 


280 


,eeeeee 


ie) 


o 
=~ 


Z 


) ee noneaaa sg 
2.0 
REL.TIME <HOURS> 


1.5 


Secon 
\ 
\ 
1.0 


) Ovo 00°0 Ov 0- 08° 0- 02'I- 
(%) NIBYLS YSAy 


Figure 11.22 Development of Radial Strain - Test MC-7.13 
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Figure 11.23 Development of Radial Strain - Test MC-7.13 
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Figure 11.24 Development of Radial Strain - Test MC-7.13 
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